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Abstract—The concepts of Danckwerts about the degrees of mixing and segregation are 
extended to the case of a continuous flow system with an arbitrary but known residence time 
distribution. For this purpose a life-expectation distribution is defined in addition to the age 
distribution. Further, a condition of maximum mixedness (minimum segregation) is defined 
for such a system. This condition, and the condition of complete segregation introduced by 
DANCKWERTS, are two opposite extremes. When the system is a reactor in which a chemical 
reaction of an arbitrary order takes place, the conversion can be calculated for both cases; thus 
two limits are obtained between which the conversion must lie. 


Résumé—Les concepts de Danckwerrrs relatifs aux degrés de mélange et de ségrégation, sont 
étendus au cas d'un systéme a écoulement continu avec une répartition de temps de séjour 
arbitrairement choisie, mais connue. Dans ce but une distribution de durée de vie probable 
est définie en plus d'une distribution @age. De plus lauteur définit pour ce systéme, une con- 
dition de mélange maximum (ségrégation minimum). Cette condition est Topposé extréme 
de la ségrégation compléte introduite par DaNcKwerts. Quand le systéme est un réacteur dans 
lequel se produit une réaction chimique d’ordre arbitraire la conversion peut étre calculée pour 
les deux cas. L’auteur obtient ainsi deux limites entre lesquelles se situe la conversion. 


Zusammenfassung Die Auffassung von Danckwerrs iiber die Grade der Mischung und 
Entmischung wird auf den Fall eines kontinuierlich durchflossenen Systems mit beliebiger aber 
bekannter Verweilzeitverteilung ausgedehnt. Hierzu wird eine Verteilung der Lebenserwartung 
zusitzlich zu der Altersverteilung definiert. Weiterhin wird eine Bedingung fiir die maximale 
Mischung (minimale Entmischung) fiir ein derartiges System definiert. Diese Bedingung und die 
Bedingung vollstiindigen Entmischung, wie sie von Danckwerrs eingefiihrt wurde, sind zwei 
entgegengesetzte Extreme. Wenn das System ein Reaktor ist, in dem eine chemische Reaktion von 
willkiirlicher Ordnung stattfindet, so kann fiir beide Fille der Umsatz berechnet werden. Man 
erhilt so zwei Grenzen, zwischen denen der Umsatz leigen muss. 


1 INTRODUCTION idealized categories. Then it is necessary to 


A PROBLEM frequently encountered in the design 
of continuous reactors for homogeneous reactions 
is the computation of the conversion from 
kinetic data obtained in batch experiments. For 
a number of simple cases, e.g. an ideal piston flow 
reactor, an ideal mixer or a series of ideal mixers, 
computation is possible and straightforward. 
In many cases, however, the reactor cannot be 


considered to fall exactly into one of these 


know more about the mixing within the system. 
A great deal of information about the state of 
mixing of the system can be obtained from the 
residence time distribution, which can often be 
found experimentally. For first order reactions 
knowledge of the residence time distribution 
suffices to compute the conversion, but with 
reactions of an order different from one, the 
residence time distribution does not determine 
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the conversion unambiguously. Danckwerts [1] 
elucidates the effect of the state of mixing on the 
average rate of reaction for a second order 
reaction. He considers a “ well-stirred”’ tank 
reactor in which the elements of fluid entering 
the reactor are uniformly dispersed through its 
volume in a time much less than the average 
residence time 7. In such a reactor the distribu- 
tion of residence times is that of an ideal mixed 
system. He then assumes two limiting cireum- 
stances. 


(a) The incoming fluid is broken up into 
discrete fragments or streaks which are 
small compared to the tank and are uni- 
formly dispersed through it, whereas the 
molecules of the fragments, which have 
entered together, remain together 
indefinitely. This is described by saying 
that the fluid remains completely 
segregated. 


(b) The inflowing material is dispersed on the 
molecular scale in a time much less than 
7; the environment of any particular 
molecule does not tend to contain an 
excess of molecules which entered at the 
same time as that molecule and the mixture 
is chemically uniform, 


The reactor output can be calculated for both 
cases, and the results show that the average rate 
of reaction of a second order reaction is greater 
when the fluid remains segregated. 

In general, mixing on the molecular scale 
between portions of reactant mixture which have 
been in the reactor for different lengths of time 
will render the output of the reactor less than if 
the material were completely segregated when 
the order of reaction is greater than unity, and 
vice versa. 

To describe the degree of mixing in a fluid, 
Danckwenrts uses the concepts of * concentration 
at a point” and “ age of the fluid at a point.” 
These imply the concentration or the age averaged 
over a region small compared to the whole system, 
and even small compared to the scale of segrega- 
tion, but large enough to contain many molecules. 


If the age « of a molecule is defined as the time 


which has elapsed since the molecule entered the 
system, it is possible to calculate the variance 
of the ages of all the molecules in the system : 
var a = (x — a) 

in which & is the mean age of all molecules which 
are at some particular moment in the system. 
The variance of the ages is the same as the mean 
square deviation of the ages of the molecules from 
the mean age. 

When for each “ point,” (defined as outlined 
above) the mean age of the molecules at that 
point is given by ap, the variance of this quantity 
can be used as a measure of the degree of segre- 
gation. This variance is: 


var ap = (ap — a)* 


in which the upper bar indicates an average over 
all * points.” In the former formula for var «, 
the upper bar indicates an average over all 
molecules. 

Referring to the two cases outlined above, 
when in a well-stirred system the mixing is uniform 
on the molecular scale, the mean age of all 
molecules at each point is the same (and equal 
to the mean age of all molecules in the system) 
and the variance between the ages of all points 
is zero. When, on the other hand, the system is 
completely segregated, the variance in ages 
between the points in the system equals the 
variance in ages between the molecules. 

Therefore, Danckwerts defines a quantity : 


J = var ap/var a = (ap — a), (x — a)? 


which he calls the degree of segregation. For a 
system which has a residence time distribution 
which equals that of an ideal mixer, the value of 
J may lie between zero (mixing on the molecular 
scale) and one (completely segregated case). 

For systems with some other residence time 
distribution the upper limit of J is again one, for 
the concept of complete segregation is then still 
applicable, but a value J = 0 is impossible 
because there must be a difference in ages at 
various points in the system. In this paper a 
condition of maximum mixedness, compatible 
with a given residence time distribution will be 
defined, which renders possible the computation 


of a lower limit of J, and of the minimum conver- 
sion for a reaction of an order higher than one 
(and equally the maximum conversion for an 
order of reaction smaller than one). 


2. Tur SIMULTANEITY OF THE MIXING 
AND THE CHEMICAL REACTION 


In the first paragraph it has been pointed out 
that two continuous reactors which have identical 
residence time distributions, can have different 
degrees of chemical conversion, when the state 
of mixing in these two reactors is different. 
Another view on the problem is explained in this 
paragraph with the aid of a simple example. 
Consider two different reactors, both consisting 
of a long narrow tube and an ideally mixed vessel 
in series. The only difference between the two 
reactors is the order of these two parts: in reactor 
I the pipe reactor precedes the mixed vessel, in 
reactor II the order is inversed (Fig. 1). The 
common residence time distribution of both 
systems is plotted in Fig. la. 


(b) 
Fig. 1 (a, b). Pipe reactor and mixed tank reactor in series, 


f(t) 


t 
Fic. 1 (ce). 


Total residence time distribution of pipe 
reactor and tank reactor in series. 
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It is easy to show that e.g. for a second order 
reaction the conversion in these two systems will 
be different, and that the conversion will be more 
complete in reactor I. | 

When the two systems are considered only as 
continuous mixing apparatus without chemical 
reaction (e.g. for diminishing concentration or 
temperature fluctuations in a steady flow of fluid) 
their effect will be exactly equal. The difficulty 
in the study of chemical reactors is that in 
general the mixing of the elements of the fluid 
occurs simultaneously with the chemical reaction. 
Now, although the overall mixing in both reactors 
I and IL is, in a way, the same, this mixing occurs 
in reactor II in an early stage of the chemical 
reaction, whereas in reactor I the mixing is done 
at a later stage. 

According to Danckwerts, reactor I has a 
higher degree of segregation than number II 
(J; >Jy). From the considerations given above 
it can be seen that the difference may be expressed 
in another way, viz. that the mixing in reactor II 
is earlier than the mixing in No. I. Going back 
to the concept of a completely segregated system, 
it will be clear that such a system can be described 
by saying that the mixing of the fluid particles is 
as late as possible, i.e. at the exit of the reactor. 
In the following paragraphs the other extreme 
ease will be described, which naturally must be 
a system in which the fluid is mixed as early as 
possible, where in both cases the residence time 


distributions are supposed to be identical. The 


last case is called the state of maximum mixedness. 
To describe these concepts in an exact way, some 
new quantities are defined in the next two 
paragraphs. 


3. RestpENcE Time DISTRIBUTION 
AND AGE DISTRIBUTION 

Throughout this paper the residence time 
distribution of the steady flow system is supposed 
to be known. It is given by the function F (t). 
The value of this function for a specified length 
of time equals the fraction of the molecules which 
remain shorter than a time t in the system. From 
this definition it is clear that F (t) is monotonously 
non-decreasing and : 


F (0) =0 (1) 
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lim F(t) =1 (2) 
When F (@) is differentiable, the function f(t) 
is given by: 


S(t) = aF (t)/dt (3) 

The function f(t) is called the residence time 

frequency function. From these functions one 

can calculate the mean residence time + and the 

variance of the residence times, which are, by 
definition : 


‘= | we or dt (4) 


and vart — or 


vart f(t)dt (5) 


The mean residence time is equal to the 
quotient of the volume of the system V and the 
volumetric flow velocity Q, which can be proved 
by the following reasoning. 

Consider the contents of the system at time 
8-0. The molecules which are then in the 
system have entered at some time before @ — 0, 
Now, in the interval between —¢— At and —1 
a quantity QM entered. Of this quantity a 
fraction F (t) had a residence time shorter than ¢, 
and therefore has already left the system at the 
moment @ = 0. Soe, the interval i— M to 

contributes a volume F Mt to 
the contents of the system at the moment @ — 0. 
Therefore, the total volume at this moment must 
be the sum of all these increments, or 


— F Qdt or 


V=Q|[l—F(oja (6) 


As this integral must be convergent, and 
1 — F(t) is monotonously non-increasing, it is 
easily proved that: 


limt [1 — F()] = 0 


By partial integration it is found from (6) that : 
[« {1 — F(t)] 

The first quantity to the right is zero for both 

limits, and therefore : 


| [1 — F(t)jt 


we (t) and so (7) 


The condition for formula (7) to be true is 
only that the quantity entering the system, 
and the quantity leaving the system at any time 
interval Af is exactly equal to QAt. It is not a 
necessary condition that there is only one entrance 
and only one exit. 

The age « of a molecule which is in the system 
at a given moment is defined as the time which has 
elapsed since it entered. The distribution function 
of the ages, ® (x), is by definition the fraction of 
the molecules at a certain moment that have ages 
smaller than «. The age frequency function is 
denoted by which is equal to d® (x) /dz. 

The molecules which have, at time @ = 0, 
an age between « and « + Ax have necessarily 
entered in the time interval — «a — Ax to — «. 
The total quantity that entered in this time 
interval was QAz, but a fraction F(x) of these 
molecules have already left the system at @ = 0 
since this fraction has a residence time shorter 
than «. The quantity of molecules with an age 
between a and x + As is therefore : 


QAx [1 — 


As the total quantity of molecules present is V, 
the fraction of the molecules which have ages 
between x and « + Az is: 


(Q/V) [1 — F(x)] Ax 


Therefore the frequency function ¢(x) is 
found to be : 


(x) = (1/7) [1 — F (8) 
The age distribution function ® is given by: 


(9) 


P(x) = | 
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The formulas (8) and (9) show that the age 
distribution at a certain moment within the 
system is uniquely determined by the residence 
time distribution. 

The mean age of the molecules is given by: 


é= [ (2) ds 


from which it is easily found with the use of (8), 
(5) and (4) that: 


= (#2 2r) or 724+ (vart)/2r (10) 


The variance of the ages is: 
var = — &)? d(x) dr 
0 


and this leads to : 


var x = (@/8r) — 477] (11) 


4 Lire Exrecration DistRiIBUTION 


In addition to the time which a molecule has 
spent already in the system, we can consider the 
time it will spend in the system from a specified 
time @ = 0 until it will leave. This time is called 
the life expectation A of the molecule. It must be 
noted that the word expectation is not used here 
in the statistical sense. It is the exact time which 
elapses between the moment @ = 0 up to the 
moment at which the molecule which we consider 
will leave the system. When the residence time 
of a molecule is denoted by ¢, it follows that for 
each molecule : 

t=-a+A 


The distribution of the life expectations of all 
the molecules in the system at a given moment is 
given by the function ¥ (A), which is the fraction 
of the molecules which are at time @ = 0 in the 
system and will, reckoned from this moment, stay 
shorter than a time A in the system. Equally, 
the life expectation frequency function is 
(A) = d¥ da. 

The life expectation distribution is found in a 
manner which is exactly analogous to that of the 
ages. 

In the interval between the times @ = A and 


0 = A + Ad a quantity QAA leaves the system. 
A fraction 1 — F(A) of these molecules were at 
time @ = 0 already in the system, and therefore 
the quantity of molecules which had, at time 
0 — 0, a life expectation between A and A + AA 
was QAA [1 — F(A)]. The ratio of this quantity 
to the total quantity present, V, is: 


¥(A) = (1/7) [1 — FQ)] (12) 


and 


W (A) — | (2°) 


So the life expectation distribution has exactly 
the same form as the age distribution. Therefore 
the mean life expectation must be equal to the 
mean age (10): 


A = 2 (var t), 2r (13) 


and equally : var A = vara (14) 


It must be mentioned here that the distribu- 
tions of x and A on the one hand and the distribu- 
tion of the residence times on the other refer to 
different populations. For the first two the 
population consists of all the molecules which are 
in the system at an arbitrarily chosen moment. 
For the residence time distribution the population 
consists of the molecules entering the system 
in an arbitrary interval of time, or of the molecules 
leaving the system in an arbitrary time interval. 
The bar above a quantity, indicating an average 
over the population, must therefore be looked 
upon with care. As an example, for every mole- 
cule, the total residence time ¢ must be equal to 
the sum of the age and the life expectation at the 
same moment : 

t=a+A 


By averaging over all molecules which are in the 
system at a given moment, it follows that : 
t=a+A 
where ¢ is the average residence time of all 
molecules in the system. From (10) and (13) 


one finds: 
= 7+ (vart)/r 


and therefore, in general t + +r 
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To evaluate +, which is usually called the mean 
residence time, the averaging of ¢ must not be 
done over all molecules present in the system, 
but over those that leave or enter. 


5. DistTRIBUTIONS FOR A SPECIFIED 
Point P 


For the following considerations the fluid in the 
system is considered to be divided into small 
volume elements called “ points,” a concept 
introduced by Danckwerrts [1]. The molecules 
within a point must be considered ideally mixed. 
Now we consider in each point at a given moment 
the age distribution of the molecules within 
this point, given by ®,(«) and the distribution 
of the life expectations of these molecules Yp (A). 
The corresponding frequency functions are the 
derivatives of these functions and are denoted 
dp (x) and (A). The point P is only defined 
at a specified moment. It does not keep its 
identity when time goes on. It can be divided 
into parts, or mixed up with others, or molecules 
can diffuse from it and to it. Therefore the 
molecules of point P can have very different life 
expectations. 

From the definitions it follows that : 


¢ (x) = (1/V) ~ (a) Up 
points 


where vp denotes the volume of point P. For 
convenience this can be written : 


(a) = (1/¥) J (18) 


since the points are supposed to be very small 
compared to the total volume. In the same way : 


(16) 


= (1/¥) bp (A) do 


Corresponding identities are valid for the distribu- 
tion functions ®(«) and ¥ (A). 


6 Tur Or MIXING OF THE 
Continuous System 


As outlined in the foregoing sections there is a 
spread between the ages of the molecules in the 
system. This spread can be characterized by the 


variance of the ages, which can be found from the 
residence time distribution with equation (11). 

For each point P, the average age at a certain 
moment can be defined by: 


ap = | da 


(17) 


Except when the system is ideally mixed on 
the molecular scale, the average ages for the 
different points are not all equal. To describe the 
spread in ages between the points, the variance of 
ap can be used, which is: 


var ap = avy (ap — a)* dv (17a) 


In this equation again the summation over all 
points is written as an integral over the volume. 
Now, Danckwerts has introduced as a measure 
for the degree of segregation, the ratio of the 
variance of the ages between the points and the 
total variance : 
J = var ap/var « (18) 
When as a third variance is introduced the 
“ variance in ages within the points,” defined by ; 


(/V) J dv | (x (19) 


it can be proved (see appendix I) that the total 
variance (equation 11) is equal to the sum of 
the variance between points (equation 17a) and 
the variance within the points (equation (19)). 
As none of these variances can be negative, the 
value of J must lie between the limits 0 and 1. 
For a completely segregated system the variance 
within the points is equal to zero, since within 
each point the ages of the molecules are all equal, 
and therefore equal to ap (see equation 19). 
It follows that in that case the variance between 
points is equal to the total variance, and J = 1. 
But a value of J = 0 is possible only for an 
ideally mixed system. In the next section the 
problem of finding the minimum value of J for a 
given residence time distribution is dealt with. 


| 
= 


7. Tur Conpition or Maximum 
MIXEDNESS 


In this paragraph a new concept is introduced, 
which is called the state of maximum mixedness. 
The significance of this concept will become clear 
when it is shown that it is the antithesis of the 
state of complete segregation, and that for this 
case J is minimum. 

The steady flow system is said to be in a state 
of maximum mixedness when two conditions are 
fulfilled. 

1. For each point, the molecules within it, 
will leave the system at the same moment, or: 


W,,(A) = for A < Ap (20) 
1 for A> Ap 
That is, the life expectation distribution 


function of the molecules at each point is a step 
function. The corresponding frequency function 
Wp (A) is a delta function: 


dp (A) = — Ap) (21) 


2. Points with equal life expectation Ap are 
mixed, or at least they have identical age distribu- 


Pp (x) = ®P (Ap, x) (22) 


The first condition is analogous to the condition 
for complete segregation, where the molecules 
at each point have entered the system at the same 
time. These conditions together express the 
requirement that all molecules, which will leave 
at the same moment and therefore will be mixed 
at the outlet, are mixed already during all the 
time they stay in the system. So it can be seen 
that when these conditions are fulfilled, the mixing 
prescribed by the (assumed) residence time 
distribution is performed as early as possible*. 


tions : 


*Here a difficulty arises with respect to the classical 
example of the ideal mixer, as for this case condition 1 
is not fulfilled. In appendix II it is shown that the 
residence time distribution of the ideal mixer is a degenerate 
case, where condition 1 can be dropped. The other 
well-known extreme case, the plug flow reactor, does 
indeed always conform to the conditions mentioned. 
In fact, the plug flow reactor is at the same time completely 
segregated and in the state of maximum mixedness, as 
this residence time distribution leaves no freedom for 
different degrees of segregation. The range of possible 
values of J has degenerated here to one value only: 
J=1. 
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In appendix II it is proved that these two con- 
ditions are sufficient for the degree of segregation 
J being minimum. The necessary conditions for 
this minimum value of J are somewhat less 
stringent. 

To show that the state of complete segregation 
and the here defined state of maximum mixedness 
are complementary we can construct two 
imaginary reactors which possess these charac- 
teristics. 

The reactor with complete segregation is 
composed of a long tube, in which piston flow is 
supposed, with a great number of side exits, 
which are placed at small intervals (see Fig. 2). 
The flow through these side exits is controlled in 
such a way that the total residence time distribu- 
tion of the system equals the given residence time 
distribution which we are studying. The volume 
of the side tubes is supposed to be negligible. 


A 


Fic. 2. 


Plug flow reactor with side exits. A case of 


complete segregation. 


As points in this system one can consider very 
short lengths of the pipe. It is clear that the age 
distribution in each point is extremely narrow, 
and that there is no mixing in the system, 
except at the exit. So the mixing of molecules of 
different age is effected as late as possible, i.e. 
in A (see Fig. 2). The position of an arbitrary 
point of the pipe can be given in two ways, viz. 
by the age « of the molecules, or by the volume v 
of the piece of pipe between the entrance and this 
point. 


Consider a small element of the pipe between 

a and « + Ax. The volume of this element is 

equal to the total volume of all the molecules with 

an age between a and « + Ax, and therefore 
equal to V¢(«). So with (8): 
Av 1 

— F (a)] Aw 

From this it follows that the relation between 

v and « is: 


(23) 
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(1 — F (a’)] da’ (24) 


Further the quantity of fluid which must be 
drawn off in the interval between « and « + Ax 
clearly must be Qf(a) Ax. By these relations 
the imaginary reactor is completely specified. 

From this first flow system, which is completely 
segregated, the second imaginary flow system, 
which is a system with maximum mixedness, is 
simply constructed by reversal of the flow. 
Then there are many entrances, and one exit 
(Fig. 3). The life expectation at any point is 
exactly determined. It is supposed that radial 
mixing is the fluid is 
immediately dispersed over the width of the pipe. 
It is easy to see that the reversal of the flow 
leaves the residence time distribution unchanged. 
Therefore this second imaginary system fulfills 
both the conditions mentioned above, and has the 
required residence time distribution. Now here 
each molecule is mixed as soon as it enters with 
the other molecules which will leave at the same 
moment. 

We now proceed to derive the explicit form of 
the function ®(Ap, «) of (22), that is, the age 
distribution of those points which have a life 
expectation equal to Ap. 

Consider again the various points of the system 
at time @ = 0. The molecules of the points which 
have life expectations equal to a specified Ap will 
leave at time @= Ap. Together with these 
molecules, however, other ones will leave, which 
enter in the time interval 0 < @ < Ap, and have 
sufliciently short residence times. The fraction 
of the leaving molecules (at time @ = Ap) which 
were not in the system at time @ = 0 equals 
F (Ap), whereas a fraction 1 — F (Ap) of these 
molecules were contained in the points with 
expectation Ap. Of these last ones, a fraction 
1 — ®,(«) was older than « at time @ = 0, and 
therefore older than «a + Ap at time @ = Ap. 
So the fraction of the leaving molecules which 
are older than « + Ap will be [1 — F (Ap)] x 
[1 — Dp (x)]. 

By definition, this fraction must be equal to 
1 — F(a + Ap), so: 


ideal, so incoming 


[1 — F (Ap)][1 — ®p(a)] = 1 — F(a + Ap) 
Solving for ®, («), 


_ F (Ap + a) — F (Ap) 

This formula is the explicit form of (22). It 

gives the common age distribution of all points 

which have a life expectation Ap, for the case of 

maximum mixedness. When F is differentiable, 
®, is clearly also differentiable, and : 


(Ap + 
1 — F (Ap) 


Using this result, it is possible to find the value 
of J for the case of maximum mixedness. To 
this end, the value of ap, i.c. the mean age of 
point P?, and the variance of zp must be computed. 


(25) 


dp (x) = (26) 


tp = dp (x) da = dz 


By partial integration and substitution of s for 
Ap + « this is found to be equal to: 


1 
| (1 — F(s)]ds (27) 


To find the numerator of J in formula (18), 
we use use (17a) in which the integration 
over the volume (summation over all points) is 
changed into an integration over Ap. The volume 
element dv must then be replaced by Vy (Ap) d Ap, 
which is equal to the joint volume of the points 
which have an expectation between Ap and 
Ap + dAp. So: 


var ap = 


fo F (s)} ds al" x 
(Ap) d Ap (28) 


In this formula, / (Ap) and a are known from (12) 
and (10) ; therefore (28) gives the variance of the 
ages between points for an arbitrary residence time 
distribution F, in the case of maximum mixedness. 
The denominator of J in (28) is given by (11), 
and thus the minimum value of J for every 
residence time distribution can be found. The 


Ap Ap 


proof that this value is indeed the minimum of J, 
is given in appendix IT. 


8 Tue CONVERSION IN A CHEMICAL 
REACTOR WITH MAXIMUM MIXEDNESS 


The amount of conversion in a reactor of 
maximum mixedness and an arbitrary residence 
time distribution can be calculated when the 
kinetics of the chemical reaction are known. 
Suppose that the incoming fluid has a concentra- 
tion ¢y, and that the reaction velocity is given by : 


r= R(e) (29) 


where R is a known function of ¢ only. 
toy 


v 


Fic. 3. Plug flow reactor with side entrances. A case of 
maximum mixedness. 


Consider now the reactor sketched in Fig. 3. 
Going along the reactor from the left to the right, 
the life expectation of the molecules A decreases, 
and becomes zero at the exit. At the left of the 
reactor the life expectation is infinite (or equal 
to the maximum residence time). Analogous to 
(23) we have here : 


Between A and A + AA a quantity Qf(A) AA is 
fed into the reactor by the side tubes. So the 


total volume flow at a point with life expectation 
A is equal to: 


{1 — F(A) JAA (30) 


[re dd 


and therefore equal to: 
Q [1 — F(A)) 


Now consider the material balance between 
A+ and A: 


Inflowing at A + Ad: 
Q [1 — F(A + Adj] e(A + Ad) 
pro unit of time 
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Outflowing at A: 
Q [1 — F(A} (A) 


Inflowing through side tubes : 
Qf (A) eg AA 


Taken away by chemical reaction : 
reaction velocity times volume, 
equal to: 
R (c)(V/r) (1) — F(A)] AA 


From this material balance one finds : 


—f(A)e(A) + [1 — F(A)] x 
(de /dd) + f(A) — R(c) [1 — =0 
or 


— C9) (31) 

This differential equation relates ¢ with A. 
The concentration at the outlet is the value of 
e for A = 0. 

The boundary condition at the left of the 
reactor for A = «© can be found by considering 
that ¢ is bounded and positive, so that for all 
normal cases de dA = 0 for A= «. The value 
of ¢ at this point can be found directly from the 
differential equation, when the limit of 
f(A) 1 — F(A) is known. When the residence 
time distribution is known only numerically, 
say from experiments with a tracer, it will appear 
that this boundary condition will give no diffi- 
culties, for the solution is nearly independent of 
the conditions at high values of A. The reason for 
this is that the very small fraction of molecules 
which have extremely high residence times does 
not affect the amount of conversion perceptibly. 
In solving the equation numerically, it is sufficient 
to choose an estimated value of ¢ for a A which is 
three or four times the mean residence time, 
and then integrate with respect to A (decreasing 
A by small steps). 


9. SOLUTION OF THE DIFFERENTIAL 
EQUATION FOR Some Cases 


a. Ideal-mixer, arbitrary reaction kinetics 
When the residence time distribution is that 
of an ideal mixer given by : 


f(t) =(1/r)e"" and F(t) =1—e*", 
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the differential equation (31) becomes: 


de 
Ble) + (1/7) (e — ea) 


The solution of this equation is de dA = 0 or 
¢ = constant. The value of ¢ is found from: 


R (c) + (1/7) (e —¢) = 0 


which is the well-known equation for a reactor 
which is ideally mixed on the molecular scale. 


b. First order reaction, arbitrary residence time 
distribution 


When the reaction is of the first order, 
R(c) = ke 
Equation (31) then reads : 


de 


which can be verified by substitution. The 
concentration at the exit is found by substitution 
of A= 0: 


= [ro e™ dt 


This expression is the same as that found when 
the system is completely segregated. In fact, for 
a first order reaction, the conversion is only 
dependent on the residence time distribution, 
whereas further differences in the degree of mixing 
have no influence. 


c. Numerical solution for two explicit residence 
time distributions 

Equation (31) has been solved for a second order 
reaction and two different residence time distri- 
butions, for which were chosen the functions : 


Sit) = (32) 


and 


(t) = (27/2) (33) 


These functions are the residence time frequency 
functions of respectively two and three well 
mixed vessels of equal volume, in series. 

If it should be known that a given reactor has 
a residence time distribution according to fy (t) 
(see above) it would still be possible to suppose 
different degrees of mixing, e.g. complete segrega- 
tion, or maximum mixedness, or that the reactor 
consists really of two well mixed vessels, in 
which case it is again possible to distinguish 
between mixing on a molecular scale in each vessel 
or complete segregation in each vessel. For all 
these cases the conversion can be calculated for a 
second order reaction. Of course many other 
intermediate cases are possible. 

Let the reaction velocity be given by: 


R (ec) = ke* 


and let the concentration at the entrance of the 
reactor be Further, let: 


(34) 


K = ke, (where K is dimensionless) (35) 


In the case of complete segregation, each fragment 
of fluid entering the reactor remains unmixed 
with other fragments in passing through the 
reactor, and its concentration decreases with time 
according to de/d@ — — ke*, which gives for an 
element of fluid with residence time ¢: 


= Kir (36) 
The concentration at the exit ¢, is then the average 
concentration of all leaving fragments of fluid, 


which is in general : 


(37) 
0 
After substitution of (82) and (36) it can be 
derived that : 
2, 4 wep 
The function Ei (x) is the exponential integral 
defined by: 


(38) 


Ei(«) = — | (e*/s) ds 
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Tables of this function can be found in [2] and [8]. 

In the case of two real vessels in series, both 
having the residence time distribution of a well 
mixed vessel, but each being completely 
segregated, the concentration after the first 
vessel (in the connecting tube) is : 


@ 
— 2 2t/r 
| 


or: 


\ K (99) 
In the same way it is found that the concentration 
at the exit of the second vessel is given by: 


From (39) and (40) the value of ¢,/cy can be 
calculated for arbitrary values of K. 

For the third case, where the reactor consists 
of two vessels in series which are ideally mixed on 
the molecular scale, the conversion is determined 
by the equations : 


= (2/7) (¢g — ¢) | (41) 
ke? = (2/1) (ce, J 
The solution for ¢, is: 


(40) 
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K 


The fourth case is the condition of maximum 
mixedness. For the given residence time distribu- 
tion (32) the differential equation (31) reduces 
to: 


(42) 


de 4A 

di — 
Let A/r = and = y 

dy 4a 

de KY 


The boundary condition is: fora = 0, dy/dx = 0. 
From this it follows that for 2 = o. 


0=Ky7¥,+2(y, —) 


From this equation y, (for 2 = «) can be 
found. For large values of a the value of y is 
nearly constant, and therefore one can start a 
numerical integration at some value of 2 which 
is so high that f(x) is negligible small, say « = 4 
or 2 = 3. The equation can then be integrated 
in small steps, and it appears that the solution 
is very nearly independent of the starting point. 
The value of y for2 = 0 (i.e. A = 0) is the required 
value of ¢,/cy which gives the concentration at 
the exit of the reactor. 


Table 1. The conversion in a reactor which has the residence time distribution of two well-mixed vessels 
in series, for different degrees of mixing, as a function keg 7. Values of ¢,/ eg 


keg t 3 5 10 20 30 30 J 
Complete segregation 0-298 0-209 0-122 0-067 0-046 0-028 1-0000 
2 real, segregated vessels Os12 0-225 O-137 0-079 0-057 0-037 0-7143 
2 real, ideally mixed vessels 0-326 0-242 0-155 0-096 0-071 0-049 0-1429 
maximum mixedness 0-252 0-166 0-106 0-081 0-9275 


Table 2. 


The conversion in a reactor which has the residence time distribution of three well-mixed 
vessels in series, for different degrees of mixing, as a function of keg r. Values of ¢,/¢9 


3 | 5 


| 


complete segregation 0-322 0-232 
3 real, segregated vessels 0-334 0-246 
3 real, ideally mixed vessels 0-357 0-275 
maximum mixedness 0-287 


10 20 30 50 J 
0-140 0-080 0-056 0-035 1-0000 
0-154 0-093 0-067 0-045 0-6250 
0-186 0-122 0-094 0-067 0-2500 
0-196 0-182 0-104 — 0-0831 


Ta. N, Zwrererine 


This integration has been carried out for 
different values of K, viz. 5, 10, 20 and 30. The 
results, together with the solutions for the other 
cases, are given in Table 1. 

Table 1 further contains the values of c¢,/e, 
for K = 3 and K = 50, for the first three cases. 
From this table it can be seen that four different 
reactors, which have identical residence time 
distributions, show different degrees of conver- 
sion. The values of the degree of segregation, 
J, are given in the last column. The value 
J = 0-0275 is the minimum value for this residence 
time distribution. 

The same calculations have been made for the 
residence time distribution f, (tf) (33), which is 
equal to the distribution for three well-mixed 
vessels in series. The formulae for the different 
cases are completely analogous (although some- 
what more complicated) to those for f, (¢). The 
results can be found in Table 2. 


10. CONCLUSION 


The considerations brought forward in this 
paper show that a knowledge of the residence 
time distribution of a continuous flow system, 
however important, is not suflicient for the 
description of the state of mixing. 

More information is, in many cases, diflicult 
to obtain. The conversion in a continuous reactor 
can be very different for different degrees of 
mixing. Probably this degree of mixing may have 
still more influence when side reactions or chain 
reactions occur. The condition of maximum 
mixedness, as defined here, may be of value in 
this field of study, giving one extreme of the state 
of mixing for every arbitrary residence time 
distribution. 


NOTATION 


concentration 
degree of segregation 
reaction velocity constant 
key + for a second order reaction 
flow velocity through the system 
residence time 
volume of the flow system 
volume of a tube reactor from the entrance up to 
a specified point 

vp = volume of a “ point ” 

a = age of a molecule 


time 
A = life expectation of a molecule 
7 = mean residence time 
exponential integral function 
residence time distribution 
mean age of the molecules with life expectation A 
age distribution 
= life expectation distribution 


Ki (x) 
F (t) 


entrance 
exit 
for a specified point 


I 
Analysis of variance 


In this appendix it is proved that the variance in ages 
of all molecules (which are at a given moment in the system) 
is equal to the sum of the variance in ages between the 
points and the variance in ages within the points, as 
defined earlier in this paper. This theorem is valid not 
only for the ages of the molecules, but for every charac- 
teristic of the individual molecules which can be expressed 
as a number, e.g. the life expectation A or the quantity 
which will be defined in appendix I, 

When 4 (2) is the frequency function of the ages of the 
molecules which are in the system at a given moment, and 
ép (x) the frequency function of the ages of the molecules 
of point P, the variances in question are : 


var x | (x — a)? d(a) dx 


var « between points : 
a iy? dv 


var x with points : 


1 | ae | ap)? dp (x) da (1, 3) 


$ (x) dx = 


(I, 4) 


(I, 5) 


( ép (a) da — 1 (for every P) 
6 


x dx (by definition) 


« 


ap adp (a) dx (by definition) 


= 1, every a) (1,8) 


(I, 6) 


(1, 7) 


€ ( 
| 
YO 
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The last equation expresses the fact that the quantity 
of molecules with ages between «and « + Ax of the system, 
which is equal to 1'é(«) Ax, is equal to the sum of the 
quantities of such molecules in all points of the system 
which is Lup ¢dp(a) As mentioned earlier, a sum 
over all points is written for convenience as an integral 
over the volume. From (I, 8) it follows (by multiplica- 
tion by « on both sides and then integrating over «) that : 


0 


Interchanging the integrations on the right hand side of 
this equation, it is found that : 


;=(1/V) | ap dv (1, 9) 


This equation states that the mean age of all molecules 
is equal to the weighed average of the mean ages of the 
points, with the volume of each point as the weight. 

With the aid of the equations (1, 4) to (I, 9) we can 
reduce the expressions for the three variances to more 
simple forms : 

var (a? — 2a + dx 
0 
a 


0 0 


(with 1, 6 and I, 4): 


x 
= | a2 d(a)dx — 24? 4 


var | (a) dx — (1, 10) 


var « between points 


| 2aap + a*) dv 


a V) | aptde — apde 


(with TI, 9) a ap? dv — 242 + 


Var a between points = (1 ~@, (1,11) 


var « within points 


ade 2a ap + ap*) dp (a) dx 
0 


a | — 2ap 
0 


dp (x) dx + ap* | ts) | 
0 


(with I, 7 and I, 5) 


=. de | — 2a%p + 


0 


Interchanging the integrations for the first term between 
the brackets, this gives : 


var « within points 


[eu | dp (x) (dv/V) — (1 [ ap? de 
Vv Vv 
(with I, 8) 


var « within points = 


b(a)dx — (1 V) ap? de (1, 12) 
4 


From the equations (1, 10), (1, 11) and (I, 12) it can be 
seen that: 


var « = var « between points + var « within points (I, 13) 


Aprenprix II 
Theorem 


The quantity J reaches its minimum (for a given 
residence time distribution) when the conditions of 
paragraph 6 are true. These conditions are sufficient but 
hot necessary. 


Proof 

Firstly, we prove that the residence time distribution 
determines uniquely the simultaneous distribution of 
a and A of the molecules which are in the system. 

Let the simultaneous frequency function of « and A 
be denoted by g(x, A), which means that a fraction 


2 (x, A) Ax AA 


of the molecules have, at moment 0 = 0, an age between 
a and a + Ax and at the same time a life expectation 
between A and A + AA. 

The quantity of molecules which have ages between 
a and « + Aa and at the same time life expectations 
between A and A + AA at time zero is the same as the 
quantity of molecules with ages between « + @ and 
a-+ 0+ Ax and life expectations between A — @ and 
A —@ + AA at time @ (when @ < A), for these are the 
same molecules, none of which has left the system. Because 
the flow is steady, all distributions are independent of time, 
therefore : 
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g(a, A) Aw MA g(a + O, A — Ax AA (il, 1) 
So, g(x, A) only cepends on the sum of « and A 
(a, A) = (a + A) 


Now the marginal distributions of g are é() and (A), 
so: 


(1, 2) 


fee. Ay dA = (11, 3) 
o 


so with (8) and (II, 2): 
| (a4 


A) dA (1/7) [1 — (HL, 4) 


(11, 5) 


Differentiation with respect to « gives: 
and so: (a, A) + A) 


From (II, 6) it is seen that the simultaneous distribution 
of « and A is completely determined by the residence time 
distribution, and therefore does not depend on the degree 
of mixing. From this it follows that also the simultaneous 
distribution of two quantities that are functions of « and 
A is completely determined, a result that is used in the 
following argument. 

Now consider a quantity ¢ which is a function of A, 
given by: 


(11, 6) 


A) da 


[et A) da 


This is the average age of the molecules having a life 
expectation A, 

For every molecule we define the quantity ¢ as the 
the given function (II, 7) in which the argument A is taken 
equal to the life expectation of this molecule (the age of 
the molecule being left out of account). For every point 
we define fp as the average of the values of £ for all 
molecules in that point. 

The analysis of variance of the quantities 2, £ and 
a — £ can be written as: 


var « within points + var « between points > 
= vara 


var ¢ within points + var ¢ between points 
= var 


var(« — within points + var(a — 
between points = var(« — £) 


Since the molecules within a given point are ideally 
mixed, the distributions of « and A within this point must 


be independent. This can be expressed by : 
A) = dp (x) dp (A) 


Therefore also x and £ must be distributed independently 
within each point. From this it follows that: 


var (a — £) within points 
var « within points + var ¢ within points 


Combining (11, 8) and (IT, 9): 


(11, 9) 


var « between points 
var « within points 


var £) within 
var € within points 


vara 


points | > (IT, 10) 


var « var (= + var 
(= £) between points + 
var within points 


The quantity at the left is the numerator of J in equation 
(18). We must prove that this quantity is minimum. 
The first two quantities of the right in eq. (11, 10) only 
depend on the residence time distribution, and not on the 
degree of mixing. Therefore they can be considered con- 
stant. The last two terms of (11, 10) can reach the value 
zero, as will be shown below, so that the minimum value 
of var « between points will be reached if, and only if 
these two quantities are both equal to zero, 


Sufficient conditions 

The quantity var ¢ within points (4th term of IT, 10) 
is zero when all molecules within a point have the same £, 
which is necessarily true when they have equal life expec- 
tations, corresponding with the first condition for maximum 
mixedness, 

The third term on the right in (II, 10) is zero when 
ap = €p for each point. This condition can be fulfilled 
owing to our choice of the function £ (A). When each point 
contains only molecules with life expectations equal to 
the value of Ap of that point, and all points with the 
same Ap have the same age distribution, each of these 
points has the sume average age ap, which must be 
equal to the average of all molecules with this life expec- 
tation, so equal to €(A) and therefore equal to &p. 

So we have proven that the two conditions for maximum 
mixedness are suflicient for reaching the minimum of the 
numerator of J, and since the denominator of J is deter- 
mined by the residence time distribution, they are 
sufficient conditions for J being minimum. 


Necessary conditions 


From (II, 10) it is clear that the necessary conditions 
for the minimum value of J are the following. 


1. All molecules within a point have equal £, that is, 
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the frequency function of € within each point is a 
6-function. When this is true, for every point fp is 
defined as the common value of & of all molecules 
within this point. 


2. Points with equal fp have equal mean ages dp. 


The first of these conditions is identical with the condi- 
tion 1 of paragraph 7, when the inverse of the function 
(A), ie. the function A (£), is single-valued. This depends 
on the given residence time distribution and can be 
ascertained by examining a plot of : 


+ 
0 


An example in which the inverse A(£) is not single- 
valued is the case of the residence time distribution 
S(t} = (l/r) exp (— t/r), (which is the distribution of a 
well-stirred vessel). For this distribution the function 
(A) is constant, viz. : 


=r (II, 12) 


Therefore the first necessary condition (all molecules of 
a point have equal £) is always satisfied, and the second 
necessary condition becomes: all points have equal mean 
ages. In the case of a vessel which is ideally mixed on the 
molecular scale, the two necessary conditions are indeed 
satisfied, whereas the conditions of paragraph 7 are not. 

The second necessary condition for minimum J, 
mentioned above, which only requires that points with 
equal € have equal ap, is clearly somewhat less stringent 
than the condition of equal a-distributions. For the 
definition of the state of maximum mixedness the more 
stringent condition is to be preferred. 
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Abstract—Equations have been developed to describe a thermal diffusion column in which 
fine packing material has been placed to reduce the convection currents. The separation calculated 
for a given column is increased when packing is introduced, in agreement with observation. 
Economic calculations show that the cost of performing a desired separation with a packed 
column is about the same as that for a column without packing. 


Résumé— Les auteurs développent des équations pour décrire une colonne a diffusion thermique 
dans laquelle ils mettent un garnissage fin pour réduire les courants de convection. La séparation 
ealeulée pour une colonne donnée augmente avec Vintroduction du garnissage comme le montre 
lexpérimentation. Du point de vue économique le prix de revient d'une séparation désirée 
est & peu prés le méme avec une colonne a garnissage ou sans garnissage. 


Zusammenfassung —Zur Beschreibung einer Thermodiffusionssiule, in der feine Fiillkérper 
zur Verringerung der Konvektion eingebracht wurden, sind Gleichungen abgeleitet worden. 
Die berechnete Trennung fiir eine gegebene Kolonne steigt durch das Einbringen der Fillkérper 
an, was in Obereinstimmung mit der Beobachtung ist. Wirtschaftlichkeitsberechnungen zeigen, 
dass die Kosten fiir cine gewiinschte Trennung mit einer gefillten Kolonne ungefihr dieselben 


wie fiir cine leere Kolonne sind. 


Tue thermal diffusion column is a mass separation 
device which appears to have some commercial 
promise for certain types of fluid systems. The 
basic thermal diffusion effect, which has been 
known for a century, is a molecular phenomenon 
in which small concentration differences arise 
under the influence of a temperature gradient. 
The thermal diffusion column, which was devised 
by CLusius and Dicken in 1938 [1], enhances 
the basic separation by means of convection 
currents, in a manner somewhat analogous to the 
way a countercurrent extraction tower produces 
concentration differences many times those 
observed in a single-stage contact. 

The principles of the column were discerned by 
a number of persons very soon after its introduc- 
tion [2, 3, 4]. The first complete presentation of 
the theory of the “open” (as contrasted to 
packed) column was that of Furry ef al. in 1989 
[5]. The column consists essentially of two 
opposing vertical plates separated by a very 
narrow open space. One plate is heated and the 
other cooled, and the thermal diffusion effect 


causes one component of, say, a binary mixture 
to diffuse toward the hot plate. At the same time, 
the density gradient which arises because of the 
temperature gradient causes smooth laminar 
convection currents up the hot plate and down the 
cold. Because of the concentration gradient set 
up by the thermal diffusion, the convection 
currents transport one component preferentially 
toward the top and thus create large concentration 
differences between the top and the bottom of the 
column. 

The distance between the hot and cold plates 
in most thermal diffusion columns using liquids 
lies in the range 0-015-0-060 in. Since the separa- 
tion obtainable decreases rapidly with increasing 
plate spacing, careful construction is obviously 
necessary. The economically optimum plate 
spacing for liquid systems is even smaller than 
this range. However, when a column with a large 
plate spacing is filled with a packing material, 
it behaves as an open column of small plate 
spacing. Thus the packed column offers the 
advantage that it requires less precision in its 
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construction than the open column. The object 
of the work reported here was to develop the 
theory of the packed thermal diffusion column 
and to determine if the introduction of packing 
alters the economic picture of thermal diffusion. 


Previous Work 


Although there is a large body of literature on 
the open column, almost no work has been 
reported on the packed column. The earliest 
precursor of the packed column appears to be the 
indication by Brewer and Bramiry [6] that 
separations are improved by the presence of a few 
baffles fastened to one plate. The first packed 
column is reported by Denyr and Burcue [7], who 
used a short column packed with glass wool to 
fractionate polymers in solution. These authors 
noted the reason that packing improves the separa- 
tion over that attained in the same column with- 
out packing. The separation achieved in a thermal 
diffusion column depends on a balance of two rates, 
the rate of mass flux by the basic thermal diffusion 
effect in the horizontal direction, and the convec- 
tion velocity in the vertical direction. Anything 
that decreases the convection velocity, such as 
inclining the column from the vertical to reduce 
the effect of gravity [8], or baffles, or packing, 
increases the relative importance of the horizontal 
flux and thus increases the steady state separation 
in the batch column. 

The only measurements of the effects of operat- 
ing variables in a packed column that have been 
reported are those of SuLLIVAN et al. [9], who 
measured the effect of plate spacing and packing 
density on separation, using glass wool packing. 
Unfortunately, not enough information was 
obtained to test the equations developed below. 


CALCULATION OF THE TRANSPORT 
EQUATION 


The calculation proceeds in several steps. 
First, the convection velocity, v, is calculated as 
a function of x. This result is used to obtain 
the horizontal mass flux, J,, then from this the 
concentration is obtained, and finally the rate 
of vertical transport, r+. 

The fluid in the column is subjected to four 
vertical forces which influence its motion, These 


are gravity, the vertical pressure gradient, fluid 
shear and the drag of the packing. It is assumed 
that the convection velocities are vertical only, 
that the pressure gradient is vertical, that there 
is no net flow through the column, that tempera- 
ture is linear in x and that density is linear 
in temperature. We follow the notation of Furry 
et al., which puts the origin of the co-ordinate 
system at the centre of the column, The differen- 
tial equation resulting from these considerations 
is 

dy gATB 

(1) 
This same equation, without the last term, is the 
starting point for the calculations on the open 
column. The last term represents the drag of 
the packing. It contains a permeability, k, which 
may be evaluated by means of isothermal flow 
measurements, using the equation 


= k(— Ap)/ul (2) 


This differs from the more traditional definition 
of permeability in that v is the actual, not the 
superficial velocity. 

The solution of equation (1), noting that the 
velocity is zero at the plates, and assuming that 
the properties of the fluid are independent of 
temperature, is 


where y = 2/4/k and Y = w/4/k. Some samples 
of the velocity profile calculated from equation 
(3) are shown in Fig. 1. The curve for Y = 0 
is the case of infinite permeability, or no packing. 
Equation (3) has been altered by the substitution 
k = w?/Y* to give an ordinate applicable to 
both packed and open columns. With this method 
of presentation, the shape and magnitude of the 
curves are general. To obtain velocities for a 
particular case, one would multiply the values of 
the ordinate by the properties and operating 
variables indicated. The same is true of Figs. 2 
and 3. The presence of the packing reduces the 
convection velocity considerably. For this reason, 
the ordinate of the curve for Y = 20 has been 
multiplied by 25 before plotting, and that for 


(3) 


sinh y 


w sinh Y 
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Fic. 1. Vertical convection velocity v as a function of 
horizontal distance for the open column (Y = 0) and for 
two conditions of packing. The cold plate is at the left 
edge and the hot plate is at the right edge of the plot. 
The ordinates for Y = 20 have been multiplied by 25 
and those for Y = 200 by 1,000 before plotting, to make 
these two curves visible with a common ordinate scale. 


Y = 200 has been multiplied by 1000 before 
plotting so that one ordinate scale could be used. 
Actually, the maximum dimensionless velocities 
involved are 3-2 10°? for Y = 0, 1-0 
for Y = 20, and 1-2 x 10-° for Y = 200, 

The horizontal mass flux rate of one constituent 
of the mixture is related to the velocity by the 
differential mass balance 


(4) 


In writing this equation, the assumptions are 
made that diffusion in the vertical direction is 
negligible, and that the concentration at a point 
is independent of time. The expression for J, 
obtained, assuming that d¢/dz is independent of 
r, is 
_ g AT Bpke de 
2p oz 
y/keoshy ‘) 
+ — etnh ¥ 

2w sinh Y (5) 

Examples of the mass flux distribution calculated 


J, 
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from equation (5) are shown in Fig. 2 for the two 
limiting values of Y, zero for the open column and 
infinity for a densely packed column. The vertical 
concentration gradient, d¢ dz, has been eliminated 
by taking the steady state value for the batch 
column, given by equation (9) with + = 0. The 
shape and magnitude of the dimensionless flux 
curve are relatively unaffected by the presence of 
packing. 

Equation (5) is based on a mass balance. In 
addition, there is the expression for J, in terms 
of the two contributions, ordinary and thermal 
diffusion, 

from which the concentration profile may be 
calculated by equating (5) and (6) and integrating. 
It is assumed that ¢(1 — ec) is independent of x. 
ATk 
2D 
ksinh y 
sinh ¥ 


+ 4/kactnh v) (7) 


Fic. 2. Horizontal mass flux rate J, as a function of 
horizontal distance for the two limiting cases Y = 0 and 
Y = © at steady state in the batch column. 


The rather odd assumption is made in this step 
that the 7’ appearing in the denominator of the 
thermal diffusion term is independent of 2, 
while dT’ /dx is replaced by AT/2w. It can be 
shown that in the case of the open column this 


The packed thermal diffusion column 


leads to the same result as obtained by Furry 
et al. by a later, more obvious simplification. 
Further, neglect of the variation of T with a is 
much less serious than the neglect of the variation 
of viscosity and diffusivity. Examples of the 
concentration profile are shown in Fig. 3 for the 
same cases as Fig. 2. The steady state batch value 
of d¢c/dz has again been used. The shape and 
magnitude of the concentration profile are 
relatively unaffected by the presence of packing. 
The rate of vertical transport of one component, 
neglecting vertical back-diffusion, is given by 
t= Bpe | ve dx (8) 

Combining (3), (7) and (8) and integrating, one 

obtains for zero throughput (the batch case) 


<4 -22 


(9) 


where 


— 8B AT? wh Boe 
6nT 


K — K, + K, 


18 15 
15? D 
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+ 


45 
etnh Y + iy? etnh? r| (11) 


(12) 


The term K of equation (9) is shown as consisting 
of two parts, K, and K,. The previous equations 
itemized here actually lead only to K,. We follow 
the procedure of Furry et al. of adding the term 
K,, which represents back-diffusion in the vertical 
direction, after the other operations are completed. 
Although back-diffusion is negligible under 
ordinary operating conditions, it is important 
in determining the optimum conditions, 

As the value of Y approaches infinity, corres- 
ponding to a fairly impermeable packing or a 
large plate spacing, the brackets in (10) and (11) 
approach the value unity, and H and K, are given 


an 


Fic. 3. Concentration c as a function of horizontal distance 
for Y = Oand Y = © at steady state in the batch column. 


simply by the product of the properties shown. 
As Y approaches zero, corresponding to a loose 
packing or to a small plate spacing, the expressions 
can be shown to approach those obtained for the 
open column. The terms in the brackets thus 
account for the relative effect of the plates in 
imposing drag on the fluid. In a densely packed 
column, the drag of the plates is of no consequence, 
while in an open column the plates present the 
only drag on the fluid. 

The transport equation, equation (9), is the 
same as the equation derived by Furry et al. 
for the open column, except that H and K are 
given by different expressions. This equation 
has been the starting point of the development of 
a number of operating equations for the open 
column. None of these derivations involved any 
assumption about the nature of H and K, and no 
characteristic of the packed column as we picture 
it is contrary to these derivations. Consequently, 
the operating equations developed and proved 
for the open column should be equally valid for 
the packed column. However, this remains to be 
proved experimentally. 


To demonstrate the differences between the two 
types of column, consider the steady state in the 
batch column. The solution of the transport 
equation for this case is [5] 


Ing = HL/K 


K, = 2w Bp De = 
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where 
= ) 
q (, —c 


If back-diffusion is ignored and the permeability 
is assumed to be low enough for the brackets in 
equations (10) and (11) to equal unity, the 
operating equation is 


/ top bottom 


5 aye DL 
29T 


Ing 


For the open column, the corresponding expression 
63 DL 
Bw* 


The two are quite similar and differ only in the 
numerical constant, the permeability, and the 
effect of plate spacing. The ratio of steady state 
separations achieved in a given column with and 
without packing is 


Ing 


(In packed 5 
(IN open 63 k 


As an example, consider a plate spacing of jin. 
and a packing of uniform spheres of diameter 
0-O2in., for which the permeability is about 
2 « 10°%em*. The ratio of separations is 4000, 
which certainly is in accord with the observation 
that the presence of packing improves the steady 
state separation, Unfortunately, it also increases 
the time it takes to reach steady state. The ratio 
of the time to reach a given fraction of the steady 
state in the packed column to that without pack- 
ing is about the square of the ratio of separations 
at steady state, 1-6 « 107, Ina practical applica- 
tion, the advantage of the former effect tends 
to be offset by the disadvantage of the latter, 
and this is undoubtedly why the economic 
calculations outlined below show little difference 
between the two types of column. 


Orptimum OPERATING CONDITIONS 


The open thermal diffusion column has received 
the attention of industrial workers, but it has had 
little use in industrial operations, reportedly 
because it is expensive. To investigate what effect 
the introduction of packing into the column has 
on the cost of separation, we have considered the 


case in which the separation and rate of through- 


put have been specified and it is desired to 
determine the optimum combination of the 
variables w, Y, B and L. 

This problem was treated by Krasny-Ercen 
in 1940 [11] for the case of gaseous isotopes in the 
open column. He noted without demonstration 
that when the properties of liquid isotopes 
were used, the column dimensions were rather 
awkward, More recently, Powers and WiLke [8] 
performed the calculations for the case of liquids 
in the middle concentration range in the open 
column, and they too found the optimum dimen- 
sions awkward. Fortunately, they found they 
were able to move away from the optimum to get 
reasonable dimensions without increasing the 
cost too much. We followed the same optimiza- 
tion procedure. Further, we used the same system 
for an example as Powers and WiLke and found 
much the same result, namely absurd optimum 
and a small increase in 
moving away from the optimum sulfliciently to 
get reasonable dimensions. 

The cost of making a separation by thermal 
diffusion consists essentially of two parts, a fixed 
charge and an operating expense. The fixed 
charge is roughly proportional to the equipment 
cost, and thus to the area of the plates, BL. The 
operating expense is chiefly heat, and the heat 
rate is obtainable from the conductivity equation, 
A BL AT /2w. Thus the cost of the separation 
per unit time is 


dimensions cost on 


cost = a, BL + a,BL/w (12) 


where a, = a,\ AT 2. The constant a, is the 
fixed charge per unit area of the plates, and a, 
is the cost of heat per B.t.u, 

In order that the separation be made, the 
variables must satisfy the appropriate operating 
equation. For the simplifying conditions of equal 
size of product streams and product compositions 
in the range of mole or weight fractions of 
0-3-0-7, this is, for continuous operation, [10] 


Hin (1-2) 


H (18) 


The terms H and K are given by equations (10), 
(11) and (12), and are rewritten here for con- 
venience as 
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K, = AT? p¢/15 D 
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The restriction imposed by the operating 
equation is incorporated by replacing L in equation 
(12) by its value from equation (13). The partial 
derivative of cost with respect to each of the three 
remaining variables is equated to zero, and these 
three equations are solved simultaneously to yield 
the optimum values of the variables. The partial 
with respect to B leads to the equation 


BH — 2-78 0 A/H, (14) 


The equation which results from taking the partial 
of cost with respect to w and eliminating B by 
using (14) is 
Ky = 
wK, —5—"w = 6 


K, (15) 


Finally, the partial with respect to Y, after 
elimination of B, gives 
/ dK, dH, 

The term w can be eliminated between (15) and 
(16) to give an equation which contains functions 
of Y only. This equation cannot be solved 
algebraically, but the results of several numerical 
solutions are plotted in Fig. 4. The optimum 
value of Y is a function of the group G = (a, a,) 
(K,/2K,)'®. At higher values of Y than shown, 
the curve fits the equation Y — 0-633 G*"', 


100 


0° t 


Fic. 4. Optimum Y as a function of cost variables and 
physical constants. Obtained from simultaneous solution 
of (15) and (16). 


It can be seen that although the cost factors 
have a fairly direct effect on the optimum Y, 
the physical properties of the system, entering 
as they do only the 1,6 power, have very little 
to do with the optimum value of Y. Further, 
numerical calculations using equation (15) have 
shown that neither the cost factors nor the 
physical properties have much influence on the 
optimum value of w. Thus, given the cost data, 
most systems and desired separations will yield 
about the same optimum values of Y and w. The 
same is true of w for the open column. B and L, 
on the other hand, will vary from case to case, 

Ostensibly, the temperature difference AT and 
average temperature 7 are also variables and 
should be included in the optimization procedure, 
However, one finds simply that it is best to operate 
at the biggest AT’ feasible, so these variables 
have been omitted from the above analysis for 
simplicity. If the viscosity and diffusivity are 
assumed to be the only properties that vary with 
temperature, the cost of operation is inversely 
proportional to about the square of AT and T 
to a rather large power which depends on the 
system. Obviously, both of these quantities are 
enhanced by a high hot plate temperature, but as 
the cold plate temperature is reduced, the AT 
is increased and the T is decreased, indicating 


' that an optimum may exist. For aqueous systems, 


the optimum ratio of the absolute hot plate 


159/60 
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temperature to that of the cold plate is a little 
over two, which is commercially unattainable. 
Thus, a practical operation would probably use tap 
water as a cooling agent and the highest tempera- 
ture heat source available which did not decom- 
pose the material or, if liquid, vaporize it. 


NuMERICAL EXAMPLE 


The procedure to be followed in using these 
results is first to evaluate G, and then to locate 
the optimum Y on Fig. 4. This is then used in 
(15) to calculate the optimum value of w, then 
(14) is used to get the optimum B, and finally 
(13) is used to get L. 

As an example, we have taken the system used 
by Powrrs and Witke. The feed is 50 mole 
percent n-heptane in benzene fed at a rate of 
1000 barrels per day, and equal product streams 
containing 30 and 70 mole per cent n-heptane are 
produced. The actual feed rate chosen affects 
only the value of B, and has no influence on the 
optimum values of the other variables. Using 
their values of system properties and cost data, 
we obtained an optimum plate spacing of 0-0087 in. 
almost the same as that for the open column, The 
optimum L and B are also similar in the packed 
and open column, as shown in Table 1. (Note 
that the values for B and “ cost ” for the packed 
column in Table 1 should be divided by the void 
fraction of the packing). The packed column 
requires one more parameter, of course, the 
optimum value of Y, which is 4-5. 


Table 1. Summary of design calculations 


i 
| None 


Arbitrarily 


0-0072 1/32 
0-56 (in.)| 17 (ft) 
i” 
8400 10,900 
2w (in.) O-0087 1/32 

L 0-41 (in.)| 3-2 (in.) 
«B (miles) 130 55 
kiem® x 108) 6-1 21 

« cost ($ day) 8100 9300 


cost ($ /day) 


* from Powers and Wiike [8) 


The optimum plate spacing is too small to be 
practical, and it is necessary to move away from 
the optimum to get a more reasonable value. If 
the plate spacing is set at 1/32in. while the other 
variables are allowed to assume their most 
optimum values subject to this limitation, a value 
of L of 3-2 in. is obtained, This is not a feasible 
column length. This same procedure in the open 
column yielded a reasonable value of 17 ft, as 
shown in Table 1, but with the packed column 
it is necessary to move farther from the optimum 
conditions. One way of doing this is to increase 
w more, but the cost is increased less if w is held 
constant at 1/32 in. and L is increased arbitrarily, 
while B and Y are allowed to assume their most 
optimum values subject to these limitations, As 
L is increased in this procedure, B and Y decrease, 
approaching in the limit a Y of zero, signifying 
no packing at all, and the ZL and B obtained in 
the open column for a w of 1,32in. As one 
example, a feasible column length of 2-5 ft is 
shown in Table 1. 

There is one interesting characteristic of the 
permeability & obtained in all of these cases, and 
that is the fact that it is high enough to restrict 
the type of packing that may be used. A _per- 
meability of 2 « 10°*cm*® may be attained with 
uniform spheres of diameter 0-02 in., for example, 
but this size sphere is much too large for the 
1/32 in, plate spacing. Smaller spheres will of 
course have a smaller permeability, and so it is 
necessary to use a different type of packing, 
such as the glass wool used successfully by 
Desye [2] and Suttivan et al. [9). 

The advantage of putting packing in a thermal 
diffusion column has been stated to be that one 
may operate at large plate spacings and get 
separations identical to those obtained at small 
plate spacings in the open column, and thus 
avoid the necessity of careful column construc- 
tion. This is true, and for research columns is a 
decided advantage. However, the calculations 
above show that for economic operation, industrial 
packed columns will have to have about the same 
plate spacing as open columns, and therefore the 
constructional advantage is lost. Further, the 
costs are roughly equivalent. The packed column 
has a slight advantage if the void fraction is high, 


VOL 
ll 
| 
wand L 
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but this is probably not worth the extra trouble of 
dealing with the packing. It must be concluded 
that for industrial operations the packed column 
has no advantage over the open column. This 
and the other conclusions above are of course 
tentative, for the packed column equations have 
not been experimentally verified. 


NOTATION 


4, ~ fixed charges per unit area of column 
a, = operating cost factor 
- cost of heat per B.t.u. 
width of column, perpendicular to # and z 
directions 
- concentration of a component of a mixture, 
mole or weight fraction 
- centreline concentration at a particular z 
diffusion coefficient 
acceleration of gravity 
constant in transport equation, equation (9) 
= horizontal mass flux of a component of the 
mixture 
k = permeability of packing 


AT 


e 


>a 


Mew 


constant in transport equation 


= vertical height of column, in the z direction 

= friction loss in packed column in isothermal! flow 
= separation factor 

T = absolute temperature 

T = arithmetic mean of plate temperatures 


difference between plate ternmperatures 
vertical convection velocity at a point 


- half the distance between the plates 


horizontal distance co-ordinate 


= 


w/ Vk 
vertical distance co-ordinate 


= thermal diffusion constant 
= — dp/dT 


concentration difference between — product 
streams in continuous column 

void fraction of packing 

thermal conductivity of fluid 


= viscosity of fluid 


density of fluid 


= rate of throughput in one section of the con- 


tinuous column, half the total feed rate in the 
example used 


= rate of vertical transport of one component 
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Abstract—Silicotungstic acid on silica gel is a very active solid acid catalyst. 

The influence of a number of variables in the preparation of the catalyst has been studied. 
A simple method for determining its activity has been developed and used to investigate the 
effect of different pretreatments given to some silica gels as well as the influence of the per- 
centage of silicotungstic acid and of particle size. 

It was found that the silica gel with the largest average pore diameter gave the most eflicient 
catalyst. This gel had the smallest specific surface area. ‘The catalytic activity increased sub- 
stantially as the catalyst particle size was reduced, indicating that only a minor part of the 
catalyst particles is accessible to the reactants. This conclusion was confirmed by the results of 
disproportionation experiments performed with cumene, The throughput could be increased 
about 10 times — the percentage of cumene conversion remaining constant — when the particle 
size was reduced by a factor of 7. It appeared that 90 per cent of the equilibrium conversion 
could still be obtained at the extremely high L.ELS.V. of 401 hrl., using a catalyst sieve fraction 
of 0-3-0-6 mm. The diameter of a silicotungstic acid molecule and the surface area of a mono- 
molecular layer formed by 1 g of this acid as calculated from crystallographic data were found 
to be 10 A and 200 m?, respectively. 


Résumé — L’acide silicotungstique sur le gel de silice constitue un catalyseur acide hétérogéne 
trés actif. 

On a examiné Vinfluence dun nombre de variables qui jouent dans la préparation du 
catalyseur. On a mis au point une méthode simple pour déterminer son activité et pour étudier 
Veffet de différents traitements préliminaires auxquels certains gels de silice ont été soumis et 
Vinfluence du pourcentage d'acide silicotungstique et de la taille des particules. 

Il a été trouve que le gel de silice ayant le plus grand diamétre moyen des particules a donné 
le catalyseur le plus efficace. Ce gel avait la surface spécitique la plus petite. L’activité catalytique 
augmentait considérablement a mesure que les dimensions des particules du catalyseur étaient 
plus petites, ce qui indique que seulement une faible partie des particules du catalyseur est 
accessible aux réactants. Cette conclusion a été confirmée par des résultats des expériences de 
disproportionnement effec uées sur le cuméne. Le débit a pu étre augmenté de dix fois — le 
pourcentage de conversion du cuméne restant constant alors que les dimensions des particules ont 
été réduites d'un facteur de 7, Ll est apparu que 90 pour cent de la conversion d'équilibre a encore 
pu étre obtenu a une vitesse de 401. cuméne |. catalyseur hr avec une fraction du catalyseur 
obtenue sur tamis de 0-3-0-6mm. Caleulés a partir de donnes crystallographiques le diamétre 
@une molécule d'acide silicotungstiqu eet la surface d'une couche monomoléculaire formée par 
1 g de cet acide, sont de 10 A et de 200 m?® respectivement. 


Zusammenfassung —Kieselwolframsiiure auf Silicagel ist ein sehr wirksamer heterogener 
saurer Katalysator. 

Der Einfluss einer Anzahl Variablen bei der Herstellung dieses Katalysators wurde untersucht. 
Eine einfache Methode zur Bestimmung seiner Wirksamkeit wurde entwickelt, mit deren Hilfe 
sowohl der Effekt verschiedener Vorbehandlungen einiger Silicagele als der Einfluss des Pro- 
zentsatzes an Kieselwolframsiure und der Korngrésse untersucht wurde. 


Es stellte sich heraus, dass das Silicagel mit dem gréssten durchschnittlichen Porendurch- 
messer den wirksamsten Katalysator lieferte. Dieses Gel besass die kleinste spezifische Oberfliche. 
Die katalytische Wirksamkeit zeigte eine erhebliche Steigung bei abnehmender Korngrésse 
des Katalysators, was darauf hindeutete, dass nur ein sehr geringer Teil der Katalysatorteilchen 
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fiir die Reaktionsteilnehmer zugiinglich ist. Diese Folgerung fand ihre Bestitigung in den 
Ergebnissen von mit Cumol durchgefiihrten Disproportionierungsversuchen. Der Durchsatz liess 
bei prozentual gleichbleibender Cumolumsetzung auf etwa das Zehnfache steigern, 


wenn die Korngrésse um einen Faktor 7 verkleinert wurde. Es ergab sich, dass beim ausserordent- 
lich hohen Durchsatz von 401, Cumol/l, Katalysator/St der. Gleichgewichtsumsatz mit einem 
Katalysator-Korngrésse von 0,3-0,6 mm doch noch bis zu 90 per cent erreicht werden konnte. 
Aus kristallographischen Daten wurden der Durchmesser der Kieselwolframsiuremolekel und die 
Oberfliche einer monomolekularen Schicht aus 1 g dieser Siure ermittelt; die gefundenen Werte 


waren 10 Abzw. 200 m2. 


INTRODUCTION 

A Great number of investigations carried out 
in the Laboratory of Chemical Technology of 
the University of Technology at Delft have 
shown that silicotungstic acid is a very active 
catalyst in many reactions proceeding under the 
influence of acid catalysts, such as polymerization 
of propene [1] and ethene [2], hydration of 
propene and ethene [3], dehydration of alcohols 
[4], alkylation of aromatic compounds with 
olefins [5], transalkylation of polyalkylbenzenes 
with benzene [5], ete. 

Evidence has been obtained [1, 2] indicating 
that a given amount of silicotungstic acid acquires 
a greater catalytic activity when applied to a 
porous carrier, VeERSTAPPEN [1] used Surinam 
bauxite for this purpose. KLiInKENBERG [2] has 
proved that increased activity can be obtained 
with silica gel as a carrier. 

Further to these investigations we have 
attempted to provide some insight into the 
structure of the silicotungstic acid on silica gel 
catalyst. 

The work described here is far from complete, 
since the study of a solid, porous catalyst has so 
many aspects that results can only be obtained 
by selecting a limited number of variables. 

The main object of the investigation was to 
find a method of reproducibly preparing a 
silicotungstic acid on silica gel type of catalyst, 
so as to enable alkylation and transalkylation 
processes to be studied with the aid of a stan- 
dardized catalyst. 


MeEaAsuRING AcTIVITY 


1. Method 

The activity of a catalyst is measured by the 
rate at which it causes a certain chemical reaction 
to proceed. It is tested preferably in the reactions 


for which the catalyst is intended. In the present 
case these were the alkylation of benzene with 
propene, of cumene with propene and transalky- 
lation of polyisopropylbenzenes with benzene. 
It has appeared, however, that a silicotungstic 
acid on silica gel catalyst that was active in the 
polymerization of olefins has a similar effect on 
the hydration of olefins, the dehydration of 
alcohols, the alkylation of aromatic compounds 
with olefins, ete. 

The activity of a catalyst accelerating a great 
many reactions, is determined preferably in that 
reaction which combines the best indication with 
a minimum of experimental difficulties. Such 
difficulties may arise from the equipment or its 
operation, the availability of base materials 
(and the purifications involved), analysis of the 
reaction product and determination of the 
conversion, serving as a measure of the activity. 

At first we determined the activity of a catalyst 
in the polymerization of ethene as follows :— 

A steel autoclave (capacity: 250ml) was 
filled with 20g of catalyst and evacuated by 
means of a water jet air pump and weighed. 
Next, 40g of ethene was introduced under 
pressure (approx. 55 atm at room temperature). 
The rotating autoclave was brought to a tempera- 
ture of 250°C within 20 min and kept at this 
temperature for exactly 1 hr. After cooling down 
in a stream of compressed air, unconverted 
ethene was drawn off and the autoclave was 
weighed again. The ethene conversion was 
calculated from the ratio of ethene converted 
to total quantity of ethene introduced. This 
method does not take into account the loss— 
during drawing off—of any C,-fraction formed 
and of volatile C; compounds. 

A drawback of this method is that a difference 
in weight of 30-35 g is determined by weighing 
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an autoclave of 20-30 kg. The pressure drop in 
such a vessel cannot be measured with great 
accuracy either. The manometers must have a 
measuring range of at least 400 atm and they 
are subjected to temperature changes and shocks, 

Prompted by a remark made by Tame ce [6] 
and a somewhat more detailed article by Jounson 
[7], we decided, in order to circumvent the above 
difficulties, to adopt the method they propose 
for measuring the activity of cracking catalysts. 

Jounson gives the following description of the 
procedure :— The polymerization rate of propene 
is determined in an all-glass system (capacity : 
75 ml) by measuring the pressure drop with the 
aid of a mercury manometer. The reaction vessel 
is kept at 200°C, the initial propene pressure being 
193 mm Hg. The experiment is carried out with 
a 100-mg sample of cracking catalyst (14-28 mesh) 
dried at 200°C for 1 hr. 

No particulars are given about equipment and 
experimental procedure. 

The apparatus we used for the activity measure- 
ments (see Fig. 1) is made entirely from pyrex 
glass and consists of two vertical cylinders, 
placed side by side, with a capacity of about 
0-51. each; they are connected at the top by a 
short capillary (internal diameter: 3 mm), which 
can be closed by a cock. 

In the left-hand vessel, at the top, there is an 
opening through which a small cylinder of coarse 
gauze, suspended from a glass hook at the ground- 
in stopper, can be introduced into the vessel. 
The right-hand vessel is connected with a closed 
mercury manometer, about 70cm long, from 
which, by means of a vernier, the pressure in the 
apparatus can be read to within 0-1 mm. A cock 
at the top of this vessel closes the inlet to a line 
which, via a three-way cock, may be connected 
either to a vacuum pump or to a propene cylinder. 
The whole apparatus is immersed in a thermos- 
tatic oil bath up to the stopper and the two cocks. 

The activity is measured as follows :— 

The oil is brought to a temperature of 150°C, 
The cylinder, filled with the catalyst, is introduced 
into the left-hand vessel. Thereupon the apparatus 
is put under vacuum for 30min to dry the 
catalyst. After drying the apparatus is closed. 
The pressure (0-01-0-lmm of Hg) does not 


appreciably rise within a period of 15 min (the 
time generally required for one measurement). 
This is at the same time a check on possible 
leakage. 

Next the cock between the two vessels is 


closed and propene is admitted to the right-hand 


Oil 


Fic. 1. Equipment for activity measurements propene 


cylinder and vacuum pump manometer. 


vessel till the pressure has reached such a value 
that, after closing the supply line and opening 
the cock between the vessels, the pressure in the 


entire vessel, i.e. the initial polymerization 
pressure, is 200mm Hg. The pressure required 
for this purpose in the right-hand vessel has 
been accurately determined by calibration, with 
a corresponding volume of inert silica gel in the 
small cylinder. Maximum variations in the 
initial pressure were only between 199-6 and 
200-4 mm. As compared with the initial pressure 
these variations are negligible. 

As the cock between the two vessels is opened, 
a stopwatch is started and the periods during 
which the pressure in the apparatus drops, e.g. 
30-0, 40-0, etc. mm Hg below the initial pressure, 
are recorded. 

The periods required to obtain equal pressure 
drops increase very rapidly with time, because 
the catalyst surface becomes covered with 
reaction products. 


| Monometer 
| 
Pita d | VOL. 
11 
| 
Rome, 
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According to Jounson the reaction proceeds 
according to the equation :— 


da/dt =k (a —a)/(1 + br) 


where 2 = pressure drop in mm Hg ; 
t = time in sec; 
a = initial pressure in mm Hg ; 
b =a constant ; 
k = activity constant in sec™}. 


This equation is a special version of the equa- 
tion derived, among others, by Larp.er [8] for 
first order reactions occuring on the surface of a 
solid catalyst and inhibited by reaction products 
formed. and Pricnarn [9] have 
applied the same formula to the catalytic decom- 
position of N,O at a platinum surface, giving 
constant values for k and b. 

The latter are determined by plotting 2/t 
against (1,1) log [a/(a —x)], according to the 
integral formula :— 


a 1} 1 a hk 

In this way a straight line is obtained with a 
slope of 2-3 (a + 1,b), the intercept on the 
av /t-axis being — kb. The value of b is calculated 
from the slope and, next, the value of k from the 
intercept. When plotting our values of a/t 
against (1/t) log [a/(a — x)], however, we did not 
obtain a straight line. From the observations 
made, no constant values could be derived for 
k and b either graphically or by calculation. 
We therefore concluded that the proposed equa- 
tion for the reaction velocity was not applicable 
in our case. When the results obtained by 
Jounson were studied in more detail, it appeared 
that his observations, too, were not in conformity 
with the equation 


a/t = 28 (1/t) (a + 1/b) log [a/(a — —k/b. 


For, whereas according to this equation the slope 
of the line is at least 2-3 a, in all cases where 
reasonable pressure drops (e.g. down to 50 per 
cent of the initial pressure) occurred, Jounson 
obtained smaller slopes. 

Attempts to draw up an equation for the 
reaction velocity in conformity with the observa- 
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tions and from which an activity constant could 
be calculated, remained unsuccessful. We are 
therefore unable to express the activity of a 
catalyst in a simple figure. It appeared, however, 
that a catalyst causing a rapid pressure drop in 
the testing apparatus, also induced a_ high 
alkylation and dehydration velocity. 

For a number of catalysts Table 1 gives the 
periods required to reach given pressure drops. 

Before considering the experimental results, we 
shall first discuss the base materials and the 
method of preparing the catalyst. 


2. Base materials and method of preparing the catalyst 

(a) Silicotungstic acid.  Silicotungstie acid 
{10, 11] is generally considered to be represented 
by the formula H,SiW,,0,,” aq. (molecular 
weight in the presence of 5 molecules of crystalliza- 
tion water: 2969). The acid can be prepared by 
adding hydrochloric acid to a solution containing 
water glass and sodium tungstate :— 


12 Na,WO, + Na,SiO, + 26 HCl + H,SiW,,0,, 
nH,O + 26 NaCl + (11-n) H,0. 


This preparation was carried out according to 
the directions given by Boorn [10]. A repro- 
ducible product was obtained which, used as a 
catalyst under identical conditions, invariably 
produced the same result. 

After drying at 120°C for 48 hr the acid contains 
about 5 molecules of crystallization water. It is 
a stable, non-hygroscopic substance, which can 
be easily dissolved in water. 

At temperatures over 300°C, the acid begins to 
decompose into WO, and SiO,; silicotungstic acid 
itself therefore cannot be used in reactions which 
require temperatures of 300°C and higher. 

A survey of the very extensive literature on 
silicotungstic acid is given in [11]. 

(b) Silica gel. The structure of silica gel has 
been widely investigated [12]. These studies 
have revealed that silica gel is composed of small 
particles or micelles. The great porosity of silica 
gel indicates that these particles have only few 
points of contact. Thus, owing to variations in 
shape, dimensions and structure of the primary 
particles and in their packing density there are a 
great many different gels. Indeed, the term 
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Catalyst : method of preparation and remarks 


Experiment I. 

Gel 0. 

Untreated. 

20°, wt silicotungstic acid. 


Experiment II. 

Gel 0, for 4 24 hr in 4N cold 
HCl changed 4 times. 

Washed with distilled water. 

20°, wt silicotungstic acid. 


Experiment 

Gel. 0, 14 days on steam-bath 
in regularly changed 4N 
HC! until acid remained 
colourless. 

Washed with distilled water. 

20°, wt silicotungstic acid. 


Experiment IV. 

Gel 0, 4 © 24hr (inel. 4 © Shr 
boiling) in IN HSO, 
changed 4 times. 

Washed with distilled water. 

20%, wt silicotungstic acid. 

2-0-2-8 mm. 


Experiment V. 

Gel 0, same treatment as in 
exp. no. IV. 

wt silicotungstic acid. 

2-8-3-4 mm. 


Experiment VI. 

Gel 0, same treatment as in 
exp. no. IV. 

20°, wt silicotungstic acid. 

344-0 mm, 


duplicate 
t 
(sec) 


4 


Experiment VII. 
Gel 0, 4 © 24hr (incl. 4 « Shr 

boiling) in 4N HCl. 
Washed with distilled water. 
20°, wt silicotungstic acid. 
2-4-2-8 mim. 


Experiment VILL. 

Gel 0, same treatment as in 
exp. no. VIF. 

20°. wt silicotungstic acid. 

4-0 mm, 


Experiment TX. 
Gel 0, 4 © 24hr (inel. 4 Sher 
boiling in IN 
washed with distilled water. 
1g cone, H,SO, on 20g gel. 
SS sample. 


Experiment X. 

24 hr (inel. 4 «8 he 
on steam-bath) in 4N HNOs 
washed with distilled water. 

20°, wt silicotungstic acid. 


Experiment 

Gel 0, treated in a separate ex- 
periment as in exp. X, 3-2 g 
carrier gave a pressure drop 
of 25mm within 12 min. 


Experiment XIII. 


Gel 24hr in 4N cold HCI 


changed once. 
Washed with distilled water. 
20°, wt silicotungstic acid. 


Experiment XIV. 


Gel 1,8 in 4N 
changed every 2nd day. 


| Shr day on steam-bath. 


Extracted with distilled water; 
(carrier not active), 
20°, wt silicotungstic acid, 


(mm) 


60 
70 


10s 


duplicate 
t 
(sec) 


| | 
(mm)| (see) (sec) 
31 31 
30 37 | 67 67 
0 78 135 132 
155 250 
60 283 | 429 
70 496 so 785 | 760 
| so S48 | 
| | 105 
| 20 25 24 193 
60 171 169 | 1005 1959 
205 206 | 
408 408 2s 63 
| Be 150 
30 20 | 34 | 
0 | 
219 224 | 
70 391 300 | 
so 6x9 48 
| | 70 268 272 
30 25 | 25 | 80466 458 
w | - - 
60 162 166 48 48 
130 458 155 157 
70 255 260 
420 430 
30 30 | 30 Is | 9 
w 61 | él w 43 42 
60 195 wo 133 158 
70 317 S11 70 276 | wal 
80 415 80 son | 500 
| 
| 
30 33 35 38 
131 2 60 | 
60 220 217 | 270 267 
70 355 | 80 | (469 471 
80 573 
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Table 1. 


Catalyst ; method of preparation and remarks—contd. 


Duplicate 
t 
(mm) (sec) 


Experiment XV. 30 

Gel 1, 4 & 24hr in 4N cold w 
Hcl. 

Rinsed with tap water. 

Washed with distilled water. 

silicotungstic acid. 

20 2-8 mm. 


Experiment XVI. 

Gel 1, treated as in exp. XV. 
20°, wt silicotungstic acid. 
mm. 


Experiment XVII 

Gel 1, treated as in exp. XV. 
20°, wt silicotungstic acid, 
4-0 mm, 


Experiment XIX. 

Gel TL, 24hr in 4N cold HCl. 
changed once. 

Washed with distilled water. 

20°, wt silicotungstic acid. 


Experiment XX. 
Gel 11, treated as in exp. XIV. 
20°, wt silicotungstic acid. 


271 
451 


3.84g sample (containing 3-2 g 


| Experiment XXIV. 


4.809 sample (containing 3-2g 


a 
(mm) 


Experiment XXII. 

Gel I, 4 © 24 hr in 4N cold 
HCl. 

Rinsed in tap water. 

Washed with distilled water. 

2 g silicotungstic acid on 20 g 
of gel. 

3-52 g sample (containing 3-2 g 
silica gel). 


Experiment XXII. 

Gel I, treated as in exp. XXI. 

4g silicotungstic acid on 20 g 
gel. 


silica gel). 


Experiment XXIII. 

Gel I, treated as in exp. XXI. 

6 g silicotungstic acid on 20 g 
gel. 

4.16g sample (containing 3-2g 
silica gel). 


Gel L, treated as in exp. X XI. 
10 g silicotungstic acid on 20 g 45 

gel. 78 
134 
236 


silica gel). 


“silica gel” denotes at least an equally wide 
range of products as the name * hydrocarbons.” 

Consequently, if a more detailed description is 
required, some characteristics have to be 
mentioned, The most important of these in the 
case of catalysts are: the specific surface area 
S (in m?/g) and the pore volume V (in ml g), 
from which the average pore diameter can be 
calculated. For cylindrical pores d = 4 (V/S) 
10‘ A, if V and S are expressed in the terms 
indicated [13). 

Three types of medium-porous silica gel, 


denoted as gel 0, gel I and gel II, were used in 
the present investigation. 

According to information supplied by the 
manufacturer (N.V. Gembo,” Winschoten, 
Holland), the term medium-porous silica gel 


comprises a gel with the following properties : 


bulk density 
actual specific weight = 


approx. 0-50 g/cm’; 
2-20 g /cm'; 
500 m?*/g; 
75 A. 


surface area 


pore diameter 


| | Duplicate 
| | (sec) 
19 9 
30 30 | 30 65 67 
75 77 181 186 
139 43 CO 30 389 390 
246 253 «60 713 716 
427 445 | 70 
80 
30 24 24 
49 
ll 60 155 185 
70 265 267 30 22 22 
259/60 80 | 454 | 40 49 49 
- 108 106 
30 27 27 | 60 201 200 
106 106 80 O17 618 
oo Is7 — - 
70 s20 | 
80 S77 | 30 16 16 
32 32 
| 25. 30 60 60 
51 | 60 109 109 
| 92 70 | 191 193 
80 B41 336 
260 | 28 
| 80 | 446 455 — | | 
30 | 24 24 —— | 30 13 | 13 
48 48 
| 50 | 91 90 | 
60136 160 
7 269 
| so 449 
29 
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Table 2 


Specific 
surface area 
(m? of 
sample) 


Specific 
surface area 
(m? per 
quantity of 
sumple 
containing 


lg 


Pore 
volume 
(ml /g of 
sample) 


Pore volume 
(ml per 
quantity of 
sample 
containing 
Ig 
silica gel) 


Gel I 
Untreated. 


Gel 1, 24 hr in 4N cold HCI (carrier 
exp. 


Gel 4 24hr in 4N cold HCl, 

20%, wt silicotungstic acid (catalyst 
exp. XV). 

Gel I, treated as sample C. 

6g silicotugstic acid on 20g gel 

(catalyst exp. XXIII) 


Gel I, treated as sample C. 
2¢ silicontungstic acid on 20g gel 
(catalyst exp. 


Gel Il, untreated. 


Gel I, 24 hr in 4N cold HCI (carrier 
exp. XEX) 

Gel 0, untreated. 

20°, wt silicotungstic acid (catalyst 

exp. I). 


Gel 0, 4 24hr in 4N cold HCI, 
20°, wt silicotungstice acid (catalyst 
exp. 


Gel 0, 4 © 2thr (inel, 4 © Shr on 
steam bath) in 4N HNO,. 20°, wt 
silicotungstic acid (catalyst exp. X). 


silica gel) 


— 


These figures, however, merely approximate 
the average properties of the product manufac- 
tured by a given method. They do not refer to a 
particular batch. If, therefore, silica gel is used 
in such investigations as are described here, the 
relevant specific data have to be determined. 

It appeared that the silica gels we used con- 
tained components giving an alkaline reaction. 


For this reason they were given a pretreatment 
with a mineral acid, washed with water and dried 
at 120°C, 

(ce) Preparation of the catalyst. (1) Silico- 
tungstic acid is prepared according to the method 
given by Boorn. 

(2) The acid so obtained is dried for 48 hr at 
120°C, 


| Average 
7 | | sama 
| (A) 
= 
A 306 506 83 
B 417 417 0-89 | 0-89 85 
345 451 0-66 77 
1959 
D 454 OSS 77 
| | 053 55 
| | 
30 
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(3) The silicotungstic acid is dissolved in such 
an amount of water that the required quantity 
of pretreated silica gel can just be immersed in it. 
The water is evaporated in a drying oven, with 
occasional stirring of the slurry. 

The catalyst is dried for 48 hr at 120°C. 

(4) The catalyst is sieved to remove any 
particles which may have been pulverized during 
the various treatments. 


3. Results of activity measurements 


In nine measurements made on 4 g of a certain 
catalyst, although conditions were identical, 
the periods required to reach given pressure drops 
showed variations of 15 per cent. These were 
found to be due to the difference in particle size 
during sampling. In one sample small particles, 
in another large ones were predominant. After 
the sample had been divided into 3 fractions 
(2-0-2-8 mm, 2-8-3-4mm and 3-4-4-0 mm), each 
sieve fraction gave duplicate determinations in 
which variations in time averaged 2 per cent, 
with a maximum of 4 per cent. 

Table 1 gives a summary of the principal 
experimental results. Unless stated otherwise, 
in all cases 4g of the sieve fraction 2-0-2-8 mm 
was used in the activity measurements. 


4. Specific surface area and pore volume of some 
carriers and catalysts 


In order to ascertain whether the structure of 
the silica gel is modified by the acid treatment and 
how its surface area and pore volume are affected 
by the application of silicotungstic acid, we had 
to gain information concerning the specific 
surface area and the pore volume of the silica 
gel used and the catalysts prepared with it. 

Of a few samples these data were determined 
for us at Koninklijke, Shell-Laboratorium, 
Amsterdam, by means of nitrogen adsorption, 
using the B.E.T.-method. 

The results of these measurements 
summarized in Table 2. 


are 


5. Discussions and conclusions 


The data given in Table 1 permit of drawing 
the following conclusions. 


(1) It has once again been proved that the 
silica gel should be pretreated with a mineral acid 
(cf. exp. I and II). 

(2) When the silica gel is treated with hydro- 
chloric acid (exp. I), sulphuric acid (exp. IV) 
or nitric acid (exp. X), almost identical results 
are obtained when, starting from a given gel, a 
treatment with 4N acid, repeated four times, is 
applied. The effect of nitric acid seems to be 
slightly better than that of the other acids. 
Experiment XIII shows that soaking of the 
silica gel in cold 4N mineral acid for 24 hr is 
sufficient ; a prolonged acid treatment or a rise 
in temperature do not result in an increased 
activity of the catalyst. 

(8) Gel I, after being given the same treatment 
as gel 0 sample and gel II, provides the most 
active catalyst. The behaviour of gel II and gel 0 
as carriers is about the same. 

(4) Mineral acid that has remained behind in 
the gel hardly contributes to the activity of the 
catalyst. From exp. IX it appears that when 
adding 1 g of concentrated sulphuric acid to 20 g 
of gel, the resulting pressure drop is only a 
fraction of that caused by the addition of 4g 
of silicotungstic acid to 20 g of gel. 

In many experiments, among them exp. XIV, 
the carrier was found to be entirely inactive. 

(5) The smaller the particles, the greater the rate 
at which the pressure decreases. The considerable 
influence of the particle diameter indicates that 
only a small part of the catalyst particle is 
accessible to propene. 

(6) The higher the percentage of silicotungstic 
acid, the greater the activity of the catalyst. The 
catalyst with 33% of acid (cf. exp. XXIV) 
appears to possess the greatest activity. 

From the results shown in Table 2 the following 
conclusions may be drawn :— 

(1) The treatment with 4N acid does not 
appreciably influence the structure of silica gel. 
Compare e.g. experiments A and B, F and G. 

(2) Neither does the application of  silico- 
tungstic acid have an important effect. The 
measurements seem to permit of the conclusion 
that the pore diameter has slightly decreased. 
Compare experiment B with C, D and E. 
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Combining the data shown in Tables 1 and 2 
it may be concluded that the catalyst based on 
gel I owes its favourable properties to its average 
pore diameter being largest, viz. 77 A (gel 0, 
catalysts 55-59 A, gel II, 50-55 A). 

It is not the gel with the largest specific surface 
area that provides the most active catalyst, 
but that with the largest pore diameter ; for the 
specific surface area of gel I is some 100 m? 
smaller than that of gel IL. 


Continuous EXPERIMENTS 


1. Method and results 


The influence of two variables, viz: (a) the 

percentage of silicotungstic acid and (b) the 
particle size of the catalyst, has been established 
in continuous experiments, using catalysts pre- 
pared from gel I. 
(a) Influence of percentage of silicotungstic acid, 
The influence of this variable has been studied in 
the alkylation of benzene with propene. For this 
purpose a liquid mixture of pure benzene 
545°-5-50°C) and propene (98%, pure) 
was passed over a fixed catalyst bed (particle 
size 2-0-2-8 mm) in a tubular reactor. Operating 
conditions were as follows: temperature 85°C, 
liquid hourly space velocity (L.HLS.V.) 161. 
liquid feed 1. catalyst hr, benzene propene molar 
ratio 4; pressure 20 atm Abs. 


The results are given in Table 3. 

It is seen that, just as in the standard activity 
test, catalyst activity increases with silico- 
tungstic acid content. 

The experiments do not establish the extent to 

which it is possible to continue increasing the 
catalytic activity by increasing the acid concen- 
tration since under the reaction conditions applied 
total propene conversion was obtained with 
20 per cent of acid. The activity test shows that 
with 33 per cent of acid the activity was still 
increasing. This was also found by BEeneKen, 
who used the same catalyst for the dehydration 
of ethyl alcohol [4]. 
(b) Influence of particle size. The influence of 
the particle size has not been studied from the 
alkylation of benzene with propene, for the 
following reason. Working conditions should 
be such that it is still possible for the propene 
conversion to increase and decrease, which means 
low temperatures and high throughputs. The 
experiments have shown that within a very narrow 
temperature range the propene conversion 
increases from 0-100 per cent. A very accurate 
temperature adjustment is therefore necessary. 
With the highly exothermic alkylation reactions, 
however, this requirement involved so many 
difficulties that another reaction was chosen to 
study the influence of the particle size, viz. the 
disproportionation of cumene : 


Table 3. Influence of the percentage of silicotungstic acid, 


| 

Silico- Molar ratio | 

tungstic 
acid 


( % wt) 


No. of 
experiment | 


of benzene 


propene | Benzene 


| 02 
| 


0-6 


| 
| 


413 840 10-6 


404-409 730 
$86 3-93 | 73-0 | 17-5 


| 9-08 3-00 16 


Composition of reaction product and alkylate (%, wt) 


Cumene | 


18-0 68-8 


Conversion 
of propene 


Residuet 


100 - — 


20-3 Os 52 


| 


| 


59 


649 


68-0 


63 23-5 | 2-9 


os 


* di.b. = diisopropylbenzene 


+ residue = tri + tetraisopropylbenzene 


VOL. 
ll 
| 
| | 
| Polymer (%) 
| | 
| - 23 80|<o1 >98 | 
32 
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2 cumene = benzene + diisopropylbenzene 
This reaction has no appreciable thermal effect. 

Cumene was passed over a fixed catalyst bed 
in a steel tubular reactor (internal diameter 
38cm). The reaction was performed in the liquid 
phase at a temperature of 230°C, a nitrogen 
pressure of 40 atm Abs. and various throughputs. 
Two catalyst sieve fractions, viz. 0-3—-0-6 mm and 
2-8-3-4 mm were employed. 

The results given in Table 4 show that the 
disproportionation of cumene is attended by the 
formation of a small amount of triisopropyl- 
benzene according to the equations :— 

2 diisopropylbenzene =< cumene + triisopropyl- 
benzene and/or 

cumene + diisopropylbenzene = benzene + 
triisopropylbenzene. 


Besides, low throughputs give maximum 


Table 4. 


cumene conversion (46-48 per cent). This points 
to the establishment of a kind of equilibrium 
between the various isopropylbenzenes. 

As regards the influence of the particle size 
it may be concluded that, just as in the activity 
test, the use of small particles results in a greater 
conversion. With a particle size of 2-8-3-4 mm, 
a cumene conversion of about 40 per cent can 
still be reached with a L.HLS.V. of 2—41/hr1.; 
with particles of 0-3-0-6 mm, the said conversion 
can even be reached at the very high L.H.S.V. 
of 40 /hr.1. 


2. Discussions and conclusions 

When considering the catalyst particles as 
globules, with diameters of 3-1 and 0-45 mm, the 
total external surface area of a given weight of 
catalyst consisting of particles with a diameter 
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(Temperature ; 230°C ; pressure : 40 atm Abs.) 


Particle size mm 


Particle size 0-3-0-6 mm 


Composition of reaction product 
% wt) 


conversion (%,) 


Composition of reaction product 
(% wt) 


MOL. 
ll 
59/60 
= 
6 1 15:3 53-5 27-3 | 46-5 19 2 16-3 518 28-5 34 48-2 
5 2 12:7 23-8 3:3 39-8 15 16-2 51-2 29-6 30 | 48°83 
12-9 60-1 | 24-4 2-7 39-9 14 14:8 53-4 29-0 2-8 46-6 
22 <3 12-0 608 26-7 O-5 39-2 ade. 16-2 | 51-2 30-4 22 48:8 
7 13-3 59-5 25-8 20 18 15-1 53-2 20-0 2-2 46-8 
| | | 20 | 280 | 15 | 440 
12 66-6 22-1 33-4 21 58-0 26-4 21 42-0 
33 
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of 0-45 mm is 3-1/0-45 = approx. 7 times as 
large as that composed of particles 0-45 mm in 
diameter. 

Although, in view of the possible influence of 
particle size on mass transport in the catalyst 
bed, only tentative conclusions can be drawn, 
it may be stated that the reaction occurs mainly 
at the surface or in a very thin layer on the 
outside of the catalyst particles. In this connex- 
ion it is pointed out that the external surface 
area of a porous globule of silica gel is very much 
larger than that of a non-porous glass bead of 
the same diameter. 

The above shows that with the rapid reactions 
occurring under the influence of the very active 
silicotungstic acid on silica gel catalyst these 
reactions are diffusion limited. There are two 
possibilities of eliminating this inhibiting effect : 

(a) by giving the carrier the largest possible 
pores and 

(b) by giving the catalyst particles the smallest 
possible dimensions. 

In the latter case the use of a fixed catalyst 
bed imposes restrictions, because the pressure 


drop across the reactor increases according as the 
particle size decreases. For this reason it may 
be desirable to change over to a suspended 
catalyst (so-called slurry process). 
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APPENDIX 


Calculation of the diameter of a silicotungstic acid molecule 
and of the surface area that can be covered by a monomole- 
cular layer of silicotungstic acid of 1g. 


Crystalline silicotungstic acid 5 aq. possesses a body- 
centred cubic lattice. The elementary cell contains two 
acid molecules and the length of the edges is 12-11 4 0-01 A 
(sioner and Gross [14]}). The silicotungstic acid 
molecules are considered as spheres touching in the space 
lattice, then the length of the diagonal of the elementary 
cell through 0, 0, 0 and 1,1, 1 <2 = the diameter of the 
silicotungstic acid molecule = 12-1143 A; so the diameter 
of a silicotungstic acid molecule is 10-5 A. 

If the molecules are arranged as in the vertical plane 
through the diagonal mentioned above, the surface area 
of a monomolecular layer formed by 1 g of acid is :— 


(1M) N (12-11) (12-11 4/2) 1/2 AP, 


M = molecular weight of silicotungstie acid and 
N = Avogadro's number. 

By substituting the numerical values it is found that 
1 g of silicotungstic acid can cover an area of 211 m?. 

In preparing a catalyst containing 20 per cent wt of 
acid, 1g of silicotungstic acid is applied to 4g of silica 
gel; 4g of gel I have a surface area of about 1600 m?; 
this area has been determined with the aid of nitrogen 
molecules. Although silicotungstic acid molecules will be 
unable to penetrate into all places where during the 
determination of the surface area nitrogen molecules could 
be adsorbed, it seems justified to suppose that, if 20 per cent 
wt of silicotungstic acid is applied, the surface of silica gel 
is not vet fully occupied up by a monomolecular acid layer. 
This is in accordance with the observation that by applying 
a higher percentage of silicotungstic acid the activity of 
the catalyst further increased. 
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Abstract—The entrainment, wave formation and atomization of horizontal films of water, 
aqueous solutions and oils have been investigated. 

The mean film thickness measured has been compared with the theoretical film thickness 
for laminar liquid flow and a smooth gas-liquid interface, Over a wide range of conditions the mean 
actual film thickness is about 0-6 times the theoretical one. 

The critical velocities for the onset of waves and atomization were determined at different 
liquid flow rates. The critical conditions for atomization have been analysed in two ways, in a 
correlation involving the Weber number (which includes the film thickness), or in a correlation 
involving the Reynolds number for film flow (which includes the rate of liquid flow). 


Résumé—Une étude a été faite de lentrainement, de la formation d’ondes et de latomisation 
de films horizontaux d'eau, de solutions aqueuses et d’huiles. 

L’épaisseur moyenne des films mesurée a été comparée avec la valeur théorique pour le cas 
d'un écoulement laminaire de liquide et une interface lisse gaz—liquide. Dans un grand intervalle 
de conditions épaisseur moyenne réelle des films est d’environ 0,6 fois la valeur théorique. 

Les vitesses critiques pour le commencement de la formation d’ondes et de latomisation 
ont été déterminées & des vitesses différentes d’écoulement du liquide. Les conditions critiques 
datomisation ont été analysées de deux fagons: dans une corrélation comprenant le nombre de 

feber (qui implique l’épaisseur du film) ou dans une corrélation comprenant le nombre de Reynolds 
pour l'écoulement du liquide (qui implique la vitesse d’écoulement du liquide). 


Zusammenfassung—Das Mitschleppen, die Wellenbildung und die Zerstéubung horizontaler 
Hautchen von Wasser, wiissrigen Lésungen und Olen wurden untersucht. 

Die gemessene mittlere Schichtdicke wurde mit dem theoretischen Wert bei laminarer 
Strémung einer Fliissigkeit und flacher Grenzfliiche zwischen Gas und Fliissigkeit verglichen. 
Unter weit verschiedenen Bedingungen ist die mittlere wirkliche Schichtdicke etwa 0,6 der 
theoretischen. 

Die kritischen Geschwindigkeiten fiir den Beginn der Wellenbildung und Zerstiubung 
wurden bei verschiedenen Fliessgeschwindigkeiten der Fliissigkeit gemessen. Die kritischen 
Bedingungen der Zerstiubung wurden in zwei Weisen analysiert, und zwar in einer Beziehung 
mit der Weberzahl (die die Schichtdicke miteinbezieht) und in einer Beziehung mit der Reynold- 
szahl fiir Schichtstrémung (die die Strémungsgeschwindigkeit der Fliissigkeit miteinbezieht). 


1. INTRODUCTION These three stages can easily be observed in 

Wuewn a gas flows over a liquid which is initially nature on rivers or lakes (cf. Figs. 8 and 9). 
at rest or in slow motion three successive stages The case where liquid at rest or in motion is in 
of the interface can be distinguished when the contact with a flowing gas or vapour occurs in 
gas velocity is increased : various equipment in use in the field of chemical 
1. the interface is flat, engineering, e.g. in furnace tubes, transfer lines, 
2. the interface is wavy, film evaporators, condensers, packed towers and 


3. the interface is wavy and liquid is atomized gas-liquid separators. In most cases, atomization 
from the wave tops. of liquid must be prevented in order to reduce 
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Fic, 1, Wind tunnel of transparent plastic. 
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liquid entrainment as much as possible. In 
literature, theoretical and experimental data on 
wave formation are scarce and no information at 
all could be found concerning the critical velocities 
at which atomization of a gas-liquid interface 
begins. In an attempt to fill this gap in our 
knowledge a number of experiments have been 
carried out. 


Il. 

Il. 1. General 

A small wind tunnel was made of transparent 
plastic, 15 15m in cross-section and 310 em 
long, see Fig. 1. In order to avoid vibrations, the 
tunnel was set on a heavy concrete foundation, 
and the connexion leading to the fan was made 
flexible. It was carefully placed horizontally. 
The upper wall was provided with a number of 
removable covers. Air was sucked through the 
tunnel at a velocity which could be varied from 
2 to 65 m_/see by closing a by-pass valve. The air 
entered through a streamlined inlet mouth 
{(1) of Fig. 1], which had four vertical plates in 
order to free the inflowing air from large vortices. 


The flow of air was determined by means of a 
measuring orifice located between the tunnel and 
the fan. For the measurement of pressure drop, 
pressure taps were mounted in the second half 
of the tunnel with a distance of 118 em between 


them. 

Liquid was supplied through a slit (8) extending 
over the whole width of the bottom of the tunnel. 
The width of the slit was 1mm, and, in some 
experiments where liquid flow rate was high, 
12mm. It was located in the inlet mouth where 
the mean air velocity was one third of that in the 
tunnel. The liquid was entrained and accelerated 
by the incoming air until it arrived at steady 
flow conditions in the tunnel. A constant head of 
liquid was ensured by means of a constant-level 
reservoir (7) with surplus surge. The liquid flow 
rate was measured by means of a rotameter which 
was regularly calibrated by weighing the flow 
over a given time. Total liquid flow rates were 
between 2-5 kg hr and 1500kg hr. At the end 
of the tunnel the liquid was recovered through a 
perforated section (10) in the bottom of the 
tunnel ; at low liquid flow rates the liquid simply 
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drained through the perforations, but at high 
flow rates it was necessary to suck it through. 
I]. 2. Device for measurement of film thickness 

The thickness 6 of the liquid film is derived 
from the electrical resistance R of the liquid 
film between two thin platinum electrodes glued 
at some distance from each other on the bottom 
of the tunnel, see Fig. 2. By calibration of the 
device (see below) a certain resistance measured 
can be related directly to the liquid height above 
the electrodes. In order to eliminate temperature 
effects, the resistance R is compared with a 
reference resistance r of the liquid at the prevailing 
temperature. ris measured between two platinum 
plates in the draining liquid at the end of the 
tunnel {(3) in Fig. 1]. In all experiments the 
difference in temperature between the liquid at 
the measuring electrodes and that of the draining 
liquid was less than 1°C. The resistances were 
measured using a 20 million MQ meter (Electronic 
Instruments Ltd.). During experiments with 
water a Philips conductivity measuring bridge 
was used. As the specific resistance of the oils 
used was too high for reliable measurements, it 
was lowered by doping. Thus the resistance was 
reduced from about 10" to 10°-10° 2 em depend- 
ing on the viscosity of the oil. The effect of 
leakage-resistance was eliminated by a series of 
suitable measurements. 

The calibration of the measuring device was 
carried out on stagnant liquid films over the 
measuring electrodes. The measuring section of 
the tunnel was for this purpose separated from 
the rest by means of dams so that a layer of liquid 
of varying thickness could be produced by 
introducing a variable amount of liquid between 
the dams. The film thickness was then measured 
with a micrometer as well as by the resistance 
method. To avoid disturbances in the electrical 
field it was necessary to locate the dams on both 
sides at a distance of at least 8cm from the 
electrodes. It was checked that it made no 
difference to the results whether the liquid was 
stagnant or whether it was continually refreshed 
by a small overilow, thus approaching more or 
less the conditions of the actual film thickness 
measurements. 

Calibrations were performed before and after 


Thickness of electrodes 0.008 mm 
Dimensions in mm 


Fic. 2. Platinum electrodes for measuring film thickness, 


each series of measurements. With the arrange- 
ment shown in Fig. 2 the film thickness was 
found to be: 8 = constant x r/R, with the con- 
stant of the order of 0-4 for 5 up to 03m. 


For 56 >03cm the constant increased with 
increasing 6. By means of a plastic model of 
two-dimensional waves placed in the liquid at 
different heights above the electrodes it was 
checked that the real mean film thickness is 
measured, provided the wave crests are parallel 
to the direction of the electrical current. Using 


the method described, film thickness could be 
measured within +. 5 to 10 per cent, the largest 
deviations occurring at the smallest film thick- 


Ill. PROCEDURE AND 


PRELIMINARY MEASUREMENTS 


Measurements were made using water, aqueous 
solutions and oils of various viscosity. During 
each run the liquid flow rate was kept constant. 
While the air flow rate was gradually increased, 
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the air mean velocity, the liquid flow rate and the 
liquid film thickness were determined. In all our 
measurements the slip between the air and the 
liquid was so large that the liquid velocity could 
be neglected with respect to the air velocity. 

Attention was given to wave formation in a 
separate series of measurements on water and a 
mineral oil. At very low air velocities the air— 
liquid interface is flat or may show small irregular 
motions ; the air velocity at which regular two- 
dimensional waves become distinctly visible was 
taken as the critical velocity for wave formation. 
With water two sets of measurements were 
carried out, one with a stagnant film retained 
between dams at both ends of the tunnel and one 
with the film entrained by the air stream, 

At sufficiently high air velocity a situation is 
reached where atomization takes place. The onset 
of atomization could clearly be observed by 
looking through the air inlet mouth towards the 
end of the tunnel ; with the film illuminated by a 
bright flood-light, the tearing-off of drops is 
easily seen against the dark background of the 
tunnel end. By illuminating different sections 
of the tunnel, it was found that, at a certain 
velocity, atomization starts near the inlet and, 
on inereasing the velocity slightly, it spreads 
over the whole length. The velocities given in this 
paper for the onset of atomization refer to con- 
ditions where atomization was first observed 
over the whole tunnel length. 

In order to check whether entrance conditions 
have an effect on the critical velocity for 
atomization, a venturi-like constriction was built 
in the wind tunnel having its throat at the 
location of the platinum electrodes for measuring 


film thickness (Fig. 3). The critical velocities for 
atomization in this constriction were found to be 
equal to those measured in the tunnel. 

The film thickness was measured 81 em from 
the end of the tunnel and in some experiments 
also in the throat of the air inlet. At air velocities 
lower than 20m ‘sec the film thickness at the 
latter location was somewhat larger than that 
measured at the end of the tunnel. At higher 
velocities no difference was noticeable. The film 
thickness in the central part of the tunnel was 
found to be slightly larger than near the sides 
and it was observed that atomization always 
started in the central region. All 5-values were 
therefore measured in the central part of the 
tunnel. 

In the experiments with oil, the bottom of the 
tunnel was cleaned with benzene and, since the 
plastic used is well wetted by oil, a continuous 
film was obtained down to the smallest film thick- 
ness. With water, however, dry spots appeared 
and the wettability of the plastic was increased by 
coating with a 5 per cent solution of polyvinyl- 
alcohol in water and subsequent hardening of 
this coating with a 4 per cent formalin solution 
in water. 

A few additional experiments were made in 
which the bottom wall of the tunnel was 
roughened. In order to arrive at a well-defined 
roughness, plastic spheres of diameters between 
0-5 and 0-75mm were glued on the bottom, 
about 160 spheres per cm*. The relative roughness 
(« Dy with Dy — 15 em) so obtained is of the 
order of 0-004, 

The shear stress 7, at the bottom of the tunnel 
(with or without liquid film) was calculated from 
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Platinum electrodes 
fer measuring film thickness 


Fic. 3. Venturi-like constriction in wind tunnel. 
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the equilibrium of forces of that part of the tunnel 
over which the drop in pressure was measured ; 
for the three smooth walls the shear stress accord- 
ing to that measured with four smooth walls 
(Fig. 10) was taken. From 7, the coeflicient of 
friction is obtained from the formula : 


This procedure gave for the dry rough tunnel 
bottom A, = 0-028 at a Reynolds number of 
2 = 105, This value compares well with A = 0-029 
reported by Nikurapse [1] for rough pipes with a 
relative sand roughness of 0-004. 

In all experiments the total rate of liquid flow 
Q was measured. In view of a more general 
application it is desirable to know the rate of 


q8 
@ 


flow q per unit of width in the central part of the 
tunnel where the film thickness is measured. 
Since the variation of the liquid velocity distribu- 
tion over the tunnel width is not known, the 
conversion factor qB/Q has been determined 
experimentally for several of the liquids used. 
To this end, plates were set on the perforated 
section at the end of the tunnel thus permitting 
separate measurement to be made of the rate of 
flow of liquid over a width of 4-5 em in the central 
part of the tunnel where the distribution of liquid 
velocities is assumed to be constant. The result 
is shown in Fig. 4 as a function of the liquid’s 
viscosity. All the values of q reported below have 
been obtained from the measured total rate of 
flow Q using the mean conversion curve of Fig. 4. 


It was to be expected that ¢B/Q would be 


1.3 
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0.1 2 4 6 8617 2 


Fic. 4. 


4 6861 2 4 6 86100 
Viscosity, ¢P 


Conversion factor for the rate of flow per unit of width as a function of liquid viscosity. 
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greater than unity ; we have no explanation for 
the value of 0-94 obtained with the most viscous 
liquid. 

Finally, the distribution of air velocities in the 
vertical direction was determined at 270 cm from 
the tunnel inlet using a small pitot tube, in order 
to find out whether fully developed flow conditions 
are obtained. The results are shown in Fig. 5. 
For comparison the velocity distribution for 
fully turbulent flow along smooth walls according 
to von KArMAN has been added to Fig 5. The 


Distance from tunnel bottom 


agreement is reasonable. This result is in accor- 
dance with the pressure drop results reported in 
section V.2, which agree with data given in 
literature for fully developed turbulent flow in 
pipes. 


IV. Lieuips Usep In THE EXPERIMENTS 


All the relevant physical properties of the 
liquids at the operating temperatures are collected 
in the following table. The oil resistance was 
lowered (see section I1.2) by doping. 


qT 


x dry smooth tunnel bottom 


v.Karmdn's universal velocity distr 
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[+s + 5.75 log 


V=30 m/sec 
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Ke) 32 34 36 40 m/sec 
Air velocity v 


Fic. 5. Velocity distribution 270 cm from the inlet of the wind tunnel Regunne: = 3°105 (smooth walls) 
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Physical properties of the liquids used 


Operating 
temp. 


Surface 
tension 
(dyne cm) 


Mineral oil 1 

Mineral oil 2 

Mineral oil 3 

Gas oil 

Kerosine 

Water + 7 per cent butanol 
Water + 3-7 per cent butanol 
Water + 1-5 per cent butanol 
Water 


V. Discussion or THE EXPERIMENTAL 
ResuLts 
V.1. 

The critical air velocities for wave formation 
were determined with water and mineral oil 3, 
the liquid films being entrained by the air stream. 
The critical velocities for oil and water respectively, 
are shown in Figs. 6 and 7 as a function of film 
thickness. Examples of the wave patterns at 
velocities slightly above critical are shown in 
Fig. 8a and 9a. Measurements have also been 
carried out with a stagnant film of water (Fig. 7). 

The critical air velocities for atomization of the 
liquid films and for the formation of waves, 
show similar trends, as can be seen from Figs. 6 
and 7. The results obtained with mineral oil 3 
and water and with the other liquids have been 
correlated in two different ways as shown in 
Figs. 10 and 11 and are discussed below (Section 
VI. 1). 

Figs. 6 and 7 show that there is a wide gap 
between the critical air velocity for wave forma- 
tion and the critical velocity at which atomization 
starts. It is seen that both critical velocities 
tend to a constant value for very thick films 
sometimes called “deep water.” These “ deep 
water” conditions are reached at a smaller film 
thickness for atomization than for wave formation. 
Jerrreys [2] developed a theory of the formation 
of waves under “ deep water” conditions. The 
waves generated are of the gravity type and the 
effect of surface tension can then be ignored. 


Wave formation and onset of atomization 


Jerrreys found: 


Pe 

The numerical value of the constant was 
determined from observations in nature and 
Jerrreys found 47 for it. Other investigators 
report velocities that are higher than those 
following from Jerrreys’ relation ; for a survey 
of the various data see Ursex [3]. 

The critical air velocities for wave formation 
calculated from Jrerrreys’ formula have been 
included in Figs. 6 and 7 for the purpose of 
reference. Jrrrrey’s formula is not applicable 
to our experiments, the waves generated on the 
shallow films being of much shorter length, so 
that surface tension must have an effect. 

A few experiments were carried out to find 
the effect of surface-active agents on wave 
formation. With an entrained water film 0-6 mm 
thick V_, was increased from 7 m ‘sec to 10 m ‘see 
by adding 1 per cent of a synthetic wetting-agent, 
“ TEEPOL 410” (cf. Keutecan [4)). 

Measurements of wave formation on entrained 
water films by Hanrarry and Encen [5] agree 
with our results as to the order of magnitude, but 
show a wide spread. The measurement of film 
thickness in Hanrarry and ENGEN’s experiments, 
however, was not completely satisfactory, as was 
already mentioned by the authors in their paper. 

The effect of wall roughness on atomization was 
studied with mineral oil 2 at only one liquid flow 
rate of g = 0-655 cm®/sec. Both with the smooth 


Density Viscosity 
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Mean air velocity Ver 
m/sec 


x,@ Experimental for entrained film 


Atomization 


Wave formation 


4 
4 


Jeffreys’ theory 2 4 m/sec 


4 6 


Fic. 


and with the rough tunnel atomization set in 
at V = 25 m/sec. In this case the film thickness 
above the plastic beads was 0-3 mm whereas for 
the smooth bottom the film thickness was 
0-9 mm. 

No measurements on wave formation were 
carried out with the rough tunnel bottom. 


4 6 8 0mm 
Film thickness 6 


4 6810 2 


Wave formation and onset of atomization for films of mineral oil 3; y = 16 cP. 


V.2. Pressure drop measurements with smooth and 
rough tunnel wall 

Measurements on pressure drop with and with- 
out liquid film were carried out with a smooth 
tunnel bottom and with the bottom covered with 
plastic beads as described in section III. The 
results obtained with the dry, smooth tunnel 
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Film thickness 6 


Fic. 7. Wave formation and onset of atomization for water films. 


bottom have been plotted in Fig. 10 together 
with the line for smooth pipes as found in litera- 
ture (e.g. Moopy [6]) showing very good agree- 
ment. From the pressure drops measured, the 
coeflicient of friction for the tunnel bottom A, 
was calculated in the way explained in section 
III. Fig. 11a shows the values of A, obtained with 


smooth and rough bottom of the tunnel and 
various rates of flow of mineral oil 2. There was 
an appreciable spread of the results and in order 
to avoid confusion, only the mean curves are 
given. For the dry, rough tunnel bottom an 
increasing Re, GAVE aN increasing A, (Fig. 11a), 
as was found by Nikurapse [1] for this range of 
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- 10, Wall shearing stress as a function of air velocity (smooth wall). 
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Fie. 11(a) 


Reynolds number and uniform sand roughness. 
Fig. 11b shows A, as a function of ¢ for Reyvaner = 
2-5 « 10°. .This figure demonstrates clearly that 
A, is no longer affected by the roughness of the 
wall above a certain rate of liquid flow. In the 
ease considered this rate of flow is: ¢ = 0-083 
em*® sec or even less. The film thickness at this 
flow rate is 0-3. mm with smooth tunnel bottom ; 
with rough tunnel bottom and the same conditions 
the oil extended about 0-1 mm above the plastic 
beads. 


VI. ANALYSIS oF THE RESULTS ON 
ATOMIZATION AND Fiutm THICKNESS 


VI. 1. Atomization 


No theory is available for the atomization of 


liquid films. The results have therefore tentatively 
been correlated by means of two dimensionless 
parameters. It was thought that some form of 


the Weber number, expressing the ratio between 
dynamic pressure of the air stream and pressure 
set-up in the film by surface tension, would be 


of importance. 
The Weber number used was defined as: 


We = 
for the other parameter S = Vy, has been 
chosen. The correlation between these two 
parameters under the critical conditions where 
atomization is first observed is shown in Fig. 12. 
For the data pertaining to air velocities above 
25 m/sec we may represent this correlation by a 
single curve. It appears that for S > 5, We,, is 
practically constant (~ 17). For air velocities 
below 25 m/sec, We,, increases beyond the correla- 
tion. For each liquid there is an air velocity 
(V_) below which no atomization occurs, 
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mineral oil 2 ; n = 26cP 
roughness of bottom wail E 0.6mm 
Retunnel * 2.5°105 


Fie. 11(b) 


Fic. 11. Coefficient of friction at the tunnel bottom wall (smooth and rough wall) 


whatever the film thickness. As a rule of thumb: 
(m/sec) ~ } o (dyne cm). 


This is translated in Fig. 12 by vertical asymp- 
totes for each liquid. 

A difficulty in applying the above correlation 
to practical cases is that usually the film thickness 
will be unknown. A rough estimate of the film 
thickness can be made from the film thickness 
correlation presented below (Section VI. 2). 

A second analysis of the onset of atomization 
involves the tentative correlation of S = Vn,/e 
with the Reynolds number related to the liquid 
flow : 


Re, = 


The result is shown in Fig. 13 ; it has the same 
features with respect to the data at air velocities 
below 25 m/sec as the Weber correlation. For 
applications of this second correlation the flow 
q per unit of width is required and this quantity 
is more easily obtainable than the film thickness. 

The results obtained on atomization and 
pressure drop with a rough tunnel bottom (ef. 
section V.1 and V.2) suggest that the correlations 
can be applied to rough walls at least when the 
film thickness is greater than the height of the 
asperities. 

The experimental data appeared to be insuffi- 
cient to decide which of the two correlations has 
the most general significance. The correlations 
should therefore not be applied if conditions 
differ much from those of the experiments. 
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Fic. 12. Correlation for onset of atomization of liquid films. 
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Fig. 13. Correlation for onset of atomization of liquid films. 
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. 14. Correlation for the thickness of horizontal liquid films entrained by an air stream, The 
experimental data have been plotted using for 7, the shear stress at the dry and smooth tunnel wall. 
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VL2. Film thickness 


The data obtained on film thickness are 
compared with those for an ideal case, viz. the 
case of a thin film in laminar flow with a flat and 
smooth gas-liquid interface. For this case it is 
easy to derive that: 

or, written in a somewhat different form : 


do (4 Re,)'/?. 


This relationship is shown in Fig. 14 by a straight 
line. 7, is the shear stress at the (smooth) gas— 
liquid interface. 

The actual film thicknesses 5 measured have 
been correlated by means of the film-thickness 
parameter 5 (7,/p,)' */v, and the Reynolds number 
for film flow Re, in Fig. 14. 

Since the shear stress at the gas-liquid interface 
involved in the first parameter was only measured 
in a few cases, the shear stress at the dry and 
smooth tunnel wall (Fig. 10) has been taken for 
it. Fig. 14 shows that for 0-5 < Re, < 2000 all 
points are below the straight line. 

The results can also be presented as a correction 
factor on the film thickness 6, calculated for the 
ideal case mentioned. We then obtain : 


M = 0-6 for 0-5 < Re, < 2000, 
0 

Data obtained with the rough tunnel bottom 
all fall within the spread of the points of Fig. 14. 

The fact that the actual film thickness is smaller 
than 5, may be explained by assuming that the 
transport of liquid mainly takes place in the tops 
of the waves moving at relatively high speeds 
and by a difference between the shear stress at 
the dry and smooth tunnel wall taken for the 
experimental data and 75. 

In a few cases the actual values of 7, were 
measured (Fig. 11). When using these values in 
the film thickness parameter, the resulting points 
in Fig. 14 would still remain below the straight 
line calculated for the ideal case mentioned above. 
However, since no further information about the 
actual values of 7, than that given in Fig. 11 is 
available, we have based the correlation on the 
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shear stress for a smooth wall which was measured 
in our experiments (Fig. 10) and in other cases 
can be calculated from correlations found in 
handbooks. 

In the experiments with water + 1 per cent 
wetting agent the air—liquid interface was smooth 
or slightly wavy. It is seen in Fig. 14 that relevant 
points are very near to the theoretical line. For 
Re, > 1500 film flow becomes turbulent and the 
points in Fig. 14 are expected to come above the 
straight line in this region. 
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NOTATION 


B = width of tunnel (= 15 em) 
Dy = hydraulic diameter of the tunnel (= 15 em) 
g = accelerating due to gravity, taken as 
980 em /sec? 
Q — total rate of liquid flow 
q = rate of liquid flow per unit of tunnel width 
R = resistance through the liquid film over the 
measuring electrodes 
Rey = 44 4, Reynolds number for film flow 
Retunnet = V Dy/vg. Reynolds number for air flow 
r = reference electrical resistance 
S = V m/e, correlation parameter 
Vo = mean air velocity in tunnel 
V., = critical velocity for wave formation or 
onset of atomization 
V,, = critical velocity for onset of atomization for 
“deep water” 
We = Pg V2 8/0, Weber number 
y = distance from the wall 
§ = mean film thickness 
8, = film thickness for laminar flow with flat and 
smooth gas-liquid interface 
« = absolute roughness of tunnel bottom 
™m = dynamic viscosity of liquid 
A = coefficient of friction 
\ = coefficient of friction at the tunnel bottom 
(with or without liquid film) 
vg = kinematic viscosity of air 
vy, = kinematic viscosity of liquid 
Pg = density of air 
py = density of liquid 
o = surface tension 
To = Shear stress at the gas-liquid interface 
7 = Shear stress at the tunnel bottom (with or 
without liquid film) 
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Abstract—-The stirred homogeneous reactor, the solid catalytic surface and the solid fuel surface 
are shown to be governed by the same pair of algebraic equations, expressing energy and material 
conservation. A relation is derived between the temperature of the reaction region and two 
dimensionless controlling parameters, one representing the chemical loading, the other representing 
heat loss. The well-known phenomena of extinction, stable and unstable burning equilibria, 
and the existence of an upper permissible limit to the heat loss are then demonstrated. The 
paper does not present new physical insights, but aims at unifying, simplifying and clarifying 
this branch of reactor theory. 


Résume—L‘auteur montre qu'un réacteur homogéne a agitation une surface catalytique solide 
et une surface combustible solide sont régis par la méme paire d’équations algébriques exprimant 
la conservation de lénergie et la conservation de la matiére. Il en déduit une relation entre 
la température de la région réactionnelle et deux paramétres sans dimensions dont l'un représente 
la charge chimque, et lautre la perte de chaleur. Les phénoménes bien connus d’extinction, 
@équilibre stable et instable de combustion et d’existence d'un maximum de perte de chaleur 
sont ensuite démontrés. L’article n’apporte pas de nouvelles connaissances physiques mais a pour 
but @unifier, simplifier et clarifier cette partie de la théorie des réacteurs. 


Zusammenfassung —Der geriihrte homogene Reaktor, die feste Katalysatoroberfliiche und 
die feste Brennstoffoberfliiche lassen sich durch das gleiche Paar algebraischer Gleichungen 
behandeln, in denen die Erhaltung der Energie und des Stoffes zum Ausdruck kommen. 
Zwischen der Temperatur der Reaktionszone und zwei dimensionslosen Einflussgréssen, von 
denen die eine die chemische * Aufladung,” die andere den Wirmeverlust kennzeichnen, wird 
eine Beziehung abgeleitet. Die wohlbekannten Phiinomene der Ausléschung, des stabilen und 
unstabilen Brennens und der Existenz einer oberen Grenze fiir den Wiirmeverlust werden 
abgeleitet. Die Arbeit bringt keine neuen physikalischen Einsichten, sie hilft aber zur 
Vereinheitlichung, Vereinfachung und Klirung dieses Teils der Reaktortheorie. 


1. INTRODUCTION 1.2 Systems considered 


1.1 Purpose The first system is the ideal stirred reactor. 
Numerous investigations [1-11] have been made This is a vessel through which a combustible gas 
of the theory of systems involving steady exother- mixture flows steadily; mixing is supposed so 
mic chemical reaction with heat loss to the intense that the gas state (temperature and 
surroundings. The present paper concerns those composition) is uniform at all points within the 
dealt with in the first ten reference papers, Its vessel. Unlike the early treatments of [12, 13, 14, 
intention is to point out the formal similarity 15], but in accordance with those of [7, 9, 9a}, 
between the systems, to simplify and unify their heat transfer is supposed to occur from the gas 
treatment, and to present particular solutions — to the reactor walls. 

of the equations in terms of dimensionless The second system is the solid catalyst surface, 
variables. The paper is a tidying-up operation, past which reactive gas flows steadily. The surface 
and does not disclose new physical phenomena. temperature, and the composition of the gas 
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immediately adjacent to it, are supposed to have 
taken up steady values controlled by the exother- 
mic reaction at the surface and by heat loss from 
the surface ; the latter might be by radiation to 
cold vessel walls. This system has already been 
studied in [1], which also contains references to 
other early work in the field of reaction 
stability. 

The third system is the solid-fuel surface, for 
example carbon, which is suspended in a stream 
of oxidizing gas, for example air. Once again, 
heat transfer other than by convection between 
fuel and gas is allowed for. The first analysis 
of this system, but without heat loss, was by 
Wacner [16] ; analyses with heat loss have been 
made in [1, 2, 3, 4, 5, 6). 

In the last two systems the reactive region is 
the solid surface. The gas composition and 
temperature adjacent to it are assumed uniform ; 
this is often the case in practice over sufliciently 
large regions for the analysis to be valid. 

The three systems must be distinguished from 
others which exhibit similar behaviour in some 
respects but which are rendered harder to analyse 
by non-uniformity of the reaction regions. Such 
other systems are the laminar flame in a pre- 
mixed gas, in which the temperature is non- 
uniform [10, 11], and the gaseous diffusion flame 
in which the fuel-oxidant ratio is non-uniform as 
well [17, 18, 19]. 

It will be supposed below that in each case the 
reaction rate depends on gas temperature and 
composition in accordance with relatively simple 
laws, i.e. increasing according to a power law 
with temperature rise and linearly with fuel or 
oxidant concentration. The extension to other 
reaction-rate functions, for example those of the 
Arrhenius type, presents no special difficulty. 
It is in any case always possible to find an 
exponent for the power law which fits the 
experimental data over a restricted temperature 
range. 


1.3 Mathematical features 

It will appear that the mathematical problem 
reduces to the simultaneous solution of two 
algebraic equations, which may be reduced, 
for the composition-dependence assumed, to a 


single quadratic equation. Solution is therefore 
easy. Trial-and-error procedures may be necessary 
with other reaction-rate functions. 

The laminar propagating flame and the diffusion 
flame on the other hand involve the solution of 
differential equations ; it is for this reason that 
they are excluded from the present treatment. 

In each case, the condition of the reaction zone 
will be shown to depend on the values of two 
dimensionless parameters. The first of these 
involves the ratio of the mass flow rate to the 
reaction zone divided by the maximum possible 
reaction rate for the system if operating adiabati- 
cally ; it will be called the loading parameter, L. 
The second involves the ratio of the maximum 
possible heat transfer rate from the reaction zone 
to the same reference reaction rate; it will be 
called the heat loss parameter, Q. 

The solution of the problem will be seen to have 
some striking characteristics. For fixed values 
of L and Q, in general three different reaction- 
zone conditions satisfy the equations. Only 
two of these represent physically stable condi- 
tions. 

If L is too large or too small, two of the solu- 
tions become imaginary. If @ is too large, there is 
no real value of L for which three real solutions 
can be found. 


2. Turk Equations 
2.1 The stirred reactor 
It will be supposed that the volumetric reaction 
rate, g’"’ in, say, cal cm*see, can be expressed by 
the formula : 


(1) 


= Kg, 


where 


q.. is a constant of the fuel and oxidant 
streams and the reactor pressure ; 

v is the reactor volume ; 

a is the fraction of the initial fuel mass which 
is still unburned ; 

7 is the reactedness (dimensionless temperature 
rise above unburned state) 

K = (n + 1)"'' a constant which ensures 
that the maximum possible value of ¢'" is 
q,.. and n is a constant of the reaction, chosen 
so as to fit approximately the experimental data. 


Equation (1) implies that the reaction is first 
order with respect to fuel. The temperature 
dependence of the reaction rate is measured by 
the value of n. 

It is further assumed that the heat loss from the 
gas to the reactor walls in unit time, Q is given 
by 


(2) 


where m is a constant depending on the mode of 
heat transfer. Once again, the exponent is 
shown so as approximately to fit the experimental 
data over the important part of the temperature 
range. 

For constant gas specific heat, c, the steady- 
flow energy equation for the reactor can then be 
written as 


me(T, —T,)r = —Q, 


where 


(3) 


mm — mass flow rate of entering mixture 


T,, — T,, = temperature rise of mixture in 
complete adiabatic steady-flow 
reaction. 

The conservation-of-mass principle applied to 

the fuel or oxidant flowing through the reactor 
yields the equation : 


—T,) (A) 


wherein the left-hand side represents the difference 
between the inflow and the outflow, while the 
right-hand side represents the rate of disap- 
pearacne of the reactant due to chemical reaction. 


Dimensionless form of the equation. We now 
introduce the dimensionless loading and heat-loss 
parameters by the definition :— 


L=me (T, _ (5) 
and Q=Qu/Gu? (6) 


Introducing these definitions into (3) and (4), 
we obtain the pair of equations governing the 
behaviour of the stirred reactor. They are :— 


Lr=Kar" 
and L(1 — a) = Ka-" 


(7) 
(8) 
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Here we leave the matter until the other 
systems have been dealt with. It will appear 
that they yield equations which are almost 
identical with (7) and (8). 


2.2 The solid catalyst surface 


It will be supposed that the reaction rate at the 
catalyst surface, q’’, in, say, cal /em*sec, is given 
by 

= Kg, 27" (9) 
where 

Gm is the maximum possible of reaction rate at 
the catalyst surface, with the given gas stream, 
when non-convective heat loss is prevented and 
Reynolds Analogy holds (these conditions lead to 
x=1—r). gm is a property of the catalyst, 
the reactants, and the stream temperature and 
pressure. 

a is the concentration in the gas at the surface 
of the less plentiful reactant, divided by its 
concentration in the unreacted stream. 

r is the surface temperature minus the stream 
temperature divided by 7, — T,, the tempera- 
ture rise of the gas in adiabatic complete reaction. 

Reaction in the gas phase is neglected. 

Heat transfer from the catalyst to its non- 
gaseous surroundings, Q’, is supposed to obey 


the equation (10) 


where Q’,, is the heat transfer rate when the 
catalyst surface is at the temperature of the 
adiabatic combustion products. 

The catalyst surface is supposed to be engaged 
in steady heat and mass transfer with the com- 
bustible stream obeying linear laws, such that 
the convective heat transfer rate per unit area is 


gec(T, ua) 
and the reactant transfer rate per unit area is 
(1 — 2) (1) 


where 


m,, is the fuel concentration in the unreacted 
stream, 


g is the mass of stream fluid reaching tem- 
perature equilibrium with unit area of 
catalyst surface in unit time, 
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go is the mass of stream fluid reaching com- 
position equilibrium with unit area of the 
catalyst surface in unit time, and 


. according to the Chilton—Colburn analogy, 
is the Lewis No. (molecular diffusivity of 
fuel divided by thermal diffusivity of 
mixture) to the minus two-thirds power. 
When the Lewis No. equals unity, Reynolds 
Analogy holds for heat and mass transfer, 
though not necessarily for friction. 


is related to the more usual heat transfer 
coeflicient h by 


gohe (12) 


where ¢ is the specific heat of the gas at constant 
pressure. 


Conservation equations. The equations of con- 
servation of energy and material can now be 
written for an element of the catalyst surface 
in the steady state. They are respectively : 


ge(T, = (18) 


and 
Kg, %7" 
1—a)=-— 14 
| a) c(T, —T.)/ma, (14) 
wherein it is recognized that ¢(7\, — T,) is 
equal to m,, times the heat of reaction of the 
combustible, 
Dimensionless form of equations. We now in- 
troduce the dimensionless loading and heat loss 
parameters in the form 


L=ge(T, — 
= / Gm 

which, substituted in (18) and (14), lead to 
Lr=Kar"—Q7-" 


(15) 
(16) 


(17) 


-9.. Kar" (18) 


Comparison of (17) and (18) with (7) and (8) 
shows that the two pairs are identical except for 
the presence of o in (18). This difference dis- 
appears if the Lewis Number is unity, i.e. when 
the Reynolds Analogy holds. 


2.3 The solid fuel surface 


We consider a solid fuel surface suspended in 
an oxidizing gas stream ; it will be supposed, for 
simplicity, that only one oxide can be formed. 
The reaction will be be taken as first order in 
surface oxygen concentration, so that the reaction 
rate, g’', in, say, cal em*see is given by 


(19) 


where 

Gm is the maximum possible of reaction rate at 
the fuel surface with the given gas stream, when 
non-convective heat transfer is absent and 
Reynolds Analogy holds (then x 1 — 7). 
is a property of the fuel and of the gas stream 
composition, pressure and temperature ; 

21S Mos Mog. the surface fractional mass oxygen 
concentration divided by that in the gas stream ; 

is defined as (T, — (1, — T,) where 
T, is the surface temperature, 7), is the tempera- 
ture of the gas stream, and 7), is the temperature 
attained by a stoichiometric mixture of gas and 
fuel, the latter being supplied at 7). 


Equation (19) is formally identical with (9) 
Similarly, equation (10) can be taken as represent- 
ing the non-convective heat loss from the fuel 
surface, as well as that from the catalyst, while 
the convective term is again represented by, 
&e (T, T,) 

The rate of transfer of oxygen to the surface, 
m’’ per unit area, is given mass transfer theory 
(e g. [20]) as 


Mog, and Mog are as above, 
r is the stoichiometric ratio, mass of oxygen 
required for unit mass of fuel. 


This transfer rate, in the steady state, equals 
the rate of oxygen consumption, S, 


ti, = /e(T, — 


where c is the mean specific heat of the gas stream 
or products at constant pressure, and it is recog- 
nised that ¢ (7), — 7',)/mog is the heat of com- 
bustion per unit mass of oxygen. 


(21) 


VOL 
ll 
the, = Moc Mos (20) 
a1 +Mgs/T 
where 
|| 
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Conservation equations. Since (19) is identical 
with (9), equation (13) expresses the conservation- 
of-energy principle for the solid fuel surface also. 

The conservation of mass is represented by 
combining equations (20) and (21). 

After introduction of the same dimensionless 
loading and heat loss parameters from (15) and 
(16), the dimensionless energy- and mass-con- 
servation equations becomes 


Lr=Kar"—Q-” 
L _(t—a) 


o0 1+2moqgr 


(22) 


= (23) 

The former equation is identical with (7) and 
(17), which hold for the stirred reactor and for the 
catalyst surface. Equation (23) differs slightly 
both from (8) and (18), through the presence of 
the o and (1 + terms. 


2.4 Comparison of three systems 

We shall consider, for simplicity, the case in 
which the Lewis Number is unity. Then o 
disappears from equations (18) and (23). Further, 
we note that mog r is of the order of 0-1 for many 
practical systems; since, in addition, « is con- 
siderably less than unity under most conditions 
of interest, the term « mog/r will be neglected. 

These simplifications ensure that equations (7) 
and (8) now hold for each of the three physical 
systems; their solution therefore is also valid 
for the three systems. From now on, therefore, 
the stirred reactor, the catalytic surface and the 
solid fuel surface will all be simultaneously under 
discussion. 


3. SOLUTION EQUATIONS 


3.1 Analytical solution 


In practice, the operator of a combustion plant 
controls the mass flow to the reaction region, 
while the heat loss relations are fixed by the 
geometry ; the interplay of mass flow, heat flow 
and reaction then settles what temperature the 
reaction region will take up. 

Correspondingly, we shall want to use equations 
(7) and (8) to tell us what + is when L and Q 
have fixed values. It is however easier to develop 
a formula for L in terms of + and Q. This follows 
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by noting that, by eliminating K « +r” between 
(7) and (8), we find 


After substitution of (24) in (7), there results a 
quadratic equation in L, with solution : 


(24) 


+ 
+4/[1 — (2Q (1 + (1 — 7)? + 


+ (Q?/K? (1 — (25) 


3.2. Graphical representation 


Equation (25) has been evaluated for various 
values of Q. The temperature exponents used 
were: ” = 8,m = 4. The results are represented 
graphically in Fig. 1, with L as independent 
variable, + as dependent variable, and Q as 
parameter, 

It is evident that, for fixed Q, the L —+r 
relation is a closed loop; so for each value of + 
there are in general two L’s [as is deducible at 
once from the + sign in (25)]; but also, for 
each L, there are in general two 7’s (or rather 
three, since + = 0 is a solution for all L and Q). 

The area occupied by the loop diminishes as Q 
increases. At Q = 0-155 the loop encloses a tiny 
area around the point L = 0-163, + = 0-64; 
for higher values of Q there are no real pairs of 
values of L and + which satisfy the equation at 


Stable equilibrium 
above ths line 


0-054 


Q is dimensioniess heat loss parameter 


of reaction region 


Dimensionless temperature 


VOL, 
ll 
159/60 Bp 
all. 
Unstable 
(g=— 
| 
02 o4 06 08 2 
Z 
Dimensionless chemical loading 
Fig. 1. 
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For a given Q the real range of L for which 
there is a real value of + is bounded by the values 
for which the tangent to the loop is vertical. The 
corresponding points on the loops are joined by a 
broken line in Fig. 1. 

The general shape of the curves of Fig. 1 is 
typical of those cases for which n > m, i.e. for 
which reaction rate depends more strongly on 
temperature than does heat loss, as is invariably 
the case in combustion practice. 


4. Discussion oF SOLUTION 
4.1 Physical significance 

The loading parameter L, it will be recalled, 
is a dimensionless measure of the mass flow rate 
through the stirred reactor or of the mass transfer 
rate between the reacting surface and the gas 
stream. Inspection of Fig. 1 shows that, when 
L and the heat-loss parameter Q are fixed, the 
reaction region can have three possible tempera- 
tures, one of which is zero (i.e. equal to that 
of the gas stream). 

Considering for a moment the highest of these 
three temperatures, (i.e. the upper half of a loop), 
we see that over a large part of the range an 
increase of L (blowing rate) causes the tempera- 
ture to rise. This is a common phenomenon 
with burning carbon. Eventually however the 
temperature falls again, until, for an L value 
greater than that giving a vertical tangent, no 
real temperature can be found which satisfies 
the steady-burning equations: the flame has 
been extinguished. 

For fixed heat-loss parameter Q, extinction also 
occurs when L is made very small. This pheno- 
menon is also well-known: if the blowing rate 
is too low, heat losses become dominant and 
prevent continuance of steady reaction. 

If Q exceeds a critical value (0-155 in the case 
considered) no reaction is possible at all. With 
fixed reactor geometry, this condition is often 
encountered when the pressure falls; for the 
reaction-rate usually falls off more rapidly than 
the heat loss rate as the pressure falls, and so 
their ratio Q decreases. 

Q = 0 corresponds to the well-known adiabatic 
reactor studied by WacGner [16], van HeerpEn 


[12], Loncwett [14] and others. The low-L 
extinction condition has now become L = 0, and 
so has ceased to be of practical importance. 


4.2 Stability 

It will now be demonstrated that the inter- 
mediate value of + for fixed L and Q represents a 
condition which cannot be obtained in practice ; 
for, although this condition satisfies the steady- 
state equations, the equilibrium is unstable. 

Suppose that a reaction system is operating at 
such a condition, represented by U in Fig. 2. 
Now let the temperature of the system rise 
infinitesimally to U’. Will equilibrium be 


restored? The point U’ lies on a steady-state 
curve which has a higher Q than that prevailing 
in the system. The system therefore cannot 
get rid of all the heat which is developed ; its 
temperature correspondingly rises. So a small 
upward deviation from the steady-state point U 
causes the system to deviate still further upward. 


5 Stable 
Unstable 


Similar arguments show that a downward 
deviation to U” precipitates a further fall of 
temperature. The whole U branch of the loop is 
therefore unstable. 

When points on the upper (8) branch of the 
loop are considered in this manner, it is easily 
seen that a small deviation from equilibrium so 
changes the heat development rate that the 
temperature tends to be restored to its equili- 
brium value for the L and Q in question. The S 
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branch therefore represents stable equilibrium 
states. 
The trivial + = 0 solution is also stable. 


Remarks 

(a) Influence of simplifying assumptions. Equa- 
tion (25) is only valid for all three systems when 
o = 1 and mog/r = 0, and indeed when the reac- 
tion rate and heat loss depend exponentially on r. 
In most practical cases these assumptions will 
not hold exactly. It should however be obvious 
that the insertion of more realistic functions into 
the equations will not fundamentally alter the 
character of the solution. 

Extinction at the upper limit of LZ is caused by 
the reaction rate increasing more steeply than 
linearly with reactedness +; extinction at the 
lower limit of L is caused by a steeper temperature 
dependence of reaction rate than that of heat loss. 
These conditions almost invariably prevail. 

The actual run of the curves on Fig. 1, and in 
particular the upper limit of Q, naturally depend 
on the values of » and m or of other parameters 
(e.g. activation energy) describing the form of the 
reaction-rate and heat-loss functions. Each case 
has to be worked out as required. However the 
upper limit of Q is always likely to be less than 
unity. 

Realistic reaction-rate functions differ from 
the exponential ones considered in not giving zero 
values at += 0. This ensures that the Q = 0 
curve, for example, does not quite touch the line 
L = 0 but sweeps round and becomes asymptotic 
to the line + = 0. Physically, this feature permits 
spontaneous ignition at low values of L. 

Similar characteristics may appear where the 
temperature of the reservoir to which the reaction 
region loses heat is higher than that of the incom- 
ing stream. This is true of some furnace situations 
for example. Corresponding changes are easily 
made to the equations. 

(b) “Chemical resistance’ to reaction. In 
calculations of the rate of reaction at catalyst 
or solid fuel surfaces, it is common to suppose that 
mass transfer alone controls the rate, that is to say 
that the concentration of the gas-phase reactant 
is reduced to zero at the surface. We can use the 
present analysis to see whether this is true. 


_“res'stonce 
"Physical resistance” 


Total “ress? 


020 040 060 080 10 


Fie, 3. 


Figure 3 shows curves of 1 (1 — «) versus L 
for various values of Q; the values n = 8, m = 4 
have been used as for Fig. 1. Now it may easily 
be shown that 1 (1 —«) is equal to the total 
“resistance” to reaction, divided by the 
“ resistance ’ of the mass transfer process alone. 
“ Resistance” is defined as stream reactant 
concentration divided by reaction rate at the 
surface. 

Inspection of Fig. 2 shows that when heat loss 
is absent (Q = 0), 1/(1 — «) has the value unity 
at low L and increases by 15 per cent over the 
whole range of stable reaction. This 15 per cent 
is the contribution of “ chemical resistance.” 

When @ is finite, 1/(1 — x) is always greater 
than unity, and can exceed 1-5 if L is very small. 
However, over most of the stable range of burning 
1/(1 — «) does not exceed 1-25. So a calculation 
of the burning rate which entirely ignores chemical 
influences will not usually over-estimate the rate 
by more than 25 per cent. A safe general rule, 
which will often give sufficient accuracy, would be 
always to reduce the rate calculated from the 
mass-transfer-control hypothesis by 10 per cent. 

Of course different temperature dependences 
(n and m) of the chemical reaction and heat loss 
will modify the numbers in the above paragraphs. 
The larger n and m are, the smaller will be the 
amount by which 1/(1 — x) exceeds unity. If 
the chemical reaction rate does not vanish at 
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7 = 0, stable reaction is possible at much lower 
temperatures than in the present case; then 
1(1 — 2) greatly exceeds unity, and the main 
“ resistance’ to reaction is chemical. This 
matter has recently been discussed in connexion 
with ammonia synthesis by BoSxJsaKkovie and 
others [21, 22]. 


5. CONCLUSIONS 


(a) Systems involving steady flow, exothermic 
chemical reaction, and non-convective heat loss, 
take up a temperature which depends on two 
dimensionless quantities, a chemical loading and a 
heat-loss parameter. 

(b) If the rate of increase of chemical reaction 
rate with temperature itself rises with temperature, 
reaction can be extinguished by increase of the 
loading parameter. 


(c) If the rate of increase of the chemical 
reaction with temperature exceeds that of the 
non-convective heat loss, reaction can also be 
extinguished by reducing the loading. 


(d) If the heat loss parameter exceeds a 


critical value, which depends on the temperature 
dependences of chemical reaction and of heat 
loss, reaction is not possible at any value of the 
loading, unless significant chemical reaction takes 
place at the gas stream temperature. 


(e) For reactions at solid surfaces, the rate of 
reaction calculated by neglecting the * chemical 
resistance may give an over-estimate of the 
order of 25 per cent. The exact amount of the 
over-estimate depends on the reaction and heat- 
loss functions, and on the loading and heat-loss 
parameters, 
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Abstract—The mechanism of the absorption of nitrogen dioxide—nitrogen tetroxide gas mixtures 
from nitrogen into water was investigated. For this purpose a wetted wall column was used in 
which a concurrent flow of a laminar liquid film and a laminar gas core with a flat velocity profile 
could be established. These flow conditions were verified separately. 


The NO, absorption experiments were analysed on the basis of the following model : 


(a) in the gas phase, NO, and N,O, are in continuous chemical equilibrium with each other 
and they are transported from the bulk to the gas—liquid interface by molecular diffusion ; 


(b) at the interface, only N,O, is dissolved in the water ; 


(c) the diffusion of N,O, into the water is accompanied by a rapid pseudo first-order reaction 
between N,O, and water. The experimental results agreed fairly well with the proposed absorp- 
tion mechanism, and from these, data were obtained on the liquid-side kinetics of the absorption 
process. 


Résumé—On a étudié la vitesse de absorption de NO,-N,O, d'un mélange avee azote dans 
leau. Pour ce but on a utilisé une colonne & paroi mouillée, dans laquelle Peau et le gaz confluent 
en écoulement laminaire. Les conditions hydro- et aerodynamiques ont été vérifiées séparément. 


Les mesures sur l'absorption de NO, ont été analysées selon le mécanisme suivant : 


(a) dans le gaz les espéces NO, et N,O, sont transferees vers la paroi par diffusion moléculaire, 
toujours mutuellement en équilibre chimique ; 


(b) & la paroi seulement N,O, est dissolu dans l'eau ; 


(c) la diffusion du N,O, dans l'eau est accompagnée par une réaction rapide de premicr ordre 
entre le N,O, et Peau. Hl y a une correspondence satisfaisante entre le modéle proposé et les 
résultats expérimentales, dont on a pu déduire des données quantitatives sur la cinétique de 
labsorption dans le liquide. 


Zusammenfassung—Der Absorptionsmechanismus von NO,-N,O, in Wasser aus einer 
Mischung mit Stickstoff wurde untersucht. Hierzu wurde eine Kolonne mit berieselter Wand 
verwendet, in der ein laminarer Flissigkeitsfilm im Gleichstrom mit einer laminaren GasstrO6mung 
mit flachem Geschwindigkeitsprofil eingestellt werden konnte. Diese Strémungsbedingungen 
wurden getrennt verwirklicht. 


Die Versuche zur NO,-Absorption wurden mit Hilfe folgender Modellvorstellung analysiert : 


(a) in der Gasphase werden NO, und N,O, aus dem Strémungskern zur Grenzfliche durch 
molekulare Diffusion transportiert und befinden sich stindig im chemischen Gleichgewicht 
miteinander ; 


(b) an der Grenzfliche wird nur N,O, im Wasser geldst ; 


(ec) die Diffusion von N,O, in Wasser ist von einer schnelien Reaktion erster Ordnung 
zwischen N,O, und Wasser begleitet. Die Versuchsergebnisse stimmten ziemlich gut mit dem 
vorgeschlagenen Absorptionsmechanismus tiberein, Aus ihnen wurden Daten zur fliissigkeits- 
seitigen Kinetik des Absorptionsprozesses erhalten. 
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1. INTRODUCTION 


Tue absorption rate of nitrogen dioxide—nitrogen 
tetroxide gas mixtures from nitrogen into water 
plays an important role in the industrial produc- 
tion of nitric acid and the removal of NO, from 
gaseous mixtures*. 

The essential chemical reactions which accom- 
pany this absorption are : 


2NO, (N,O,) + H,O = HNO, + HNO, 
4HINO, 2NO + 2NO,(N,O,) + 2H,0 


The resulting overall reaction is represented by 
the following equation : 


3NO, (32 N,0O,) + H,O 2HNO, + NO 


The absorption of NO, in water is carried out 
in absorption towers which can be operated 
under various conditions. As a consequence, 
various designs of absorption towers are possible. 
An optimum design can only be realized, however, 
if the chemical and physical mechanisms of the 
absorption are thoroughly understood. 

During the last forty years several investigators 
have studied this process but so far their views 


on this subject are rather conflicting as appears 
from a literature survey [1]. A summary of the 
most important conclusions from the studies by 
other investigators is given below. 

CuamBers and Suerwoop [2] interpreted their 


results on the basis of the “two film theory 
and found the gas film to be rate controlling. 
However, the absorption rate was lower than 
was expected from similar absorption experiments 
with other gases. According to these authors, 
this was due to the formation of mist in the gas 
phase. In addition, they found NO to be present 
in the gas phase when a NaOH solution was 
used as the absorbing liquid ; this indicated that 
a chemical reaction between NO, and H,O 
took place in the gas phase. 

and Prince [3] and CavpLe and 

*In the gas phase NO, is in continuous equilibrium 
with N,O, according to the equation: 2NO, 
This equilibrium, which is dependent on temperature and 
pressure, is reached very rapidly. In this paper NO, will 
be referred to as the mixture of NO, and N,O,, unless 
stated otherwise. 


Densicu [4] concluded from their experiments 
that the absorption rate was controlled by the 
chemical reaction between N,O, and water in 
the liquid phase. These authors did not find an 
indication of gas phase reaction. 

Perers et al. [5] stated that the reaction 
between N,O, and H,O in the gas phase close 
to the interface was rate controlling. 

Wenpet and Picrorp [6] recently stated the 
controlling mechanism in the absorption to be the 
homogeneous liquid phase reaction of dissolved 
N,O, with water. However, their data indicate 
that the gas phase resistance also plays an 
important role in the absorption rate. 

Because of the apparent complexity of the 
problem we have especially aimed at carrying 
out the absorption experiments in an apparatus 
in which the flow conditions of both the gas and 
the liquid phase were well determined, so that a 
quantitative interpretation of the experimental 
results might be possible. In this respect, a 
favourable case would be the situation, in which 
gas and liquid, both in laminar flow, are brought 
into contact with each other for a specified period 
of time + on their way through the absorption 
unit. If, moreover, there is no velocity gradient 
perpendicular to the interface in either phase , 
it is then possible to apply a rather simple form 
of the theory of diffusion to the transfer of NO, 
from the gas to the liquid. A wetted wall column 
seemed best suited to obtain the situation 
described above, because it has three attractive 
features : 


(1) a sharp boundary between the gas and 
the liquid phase, 

(2) a well-defined absorption area, and 

(3) well-known hydrodynamic conditions of 
the falling liquid film. 


In addition, special measures were taken in 
order to insure : 

(4) cocurrent laminar flow of both gas and 
liquid, 
a flat velocity profile of the gas in the 
wetted tube, 
equal values of the surface velocity of the 
liquid film and the velocity of the gas, 
equal temperatures of gas and liquid. 


(5) 
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Before the NO, absorption experiments were 
started, a check on the flow conditions of both 
liquid and gas was desired. The laminar flow 
behaviour of the falling liquid film was verified 
by measurement of the absorption rate of pure 
CO, gas into water. The flow conditions of the 
gas were tested separately by means of absorption 
of NH, from N, into water. 


GAS MIXTURE 


Fic. 1. Sketch of absorption apparatus ; dimensions in mm. 


2. EXPERIMENTAL ParT 


2.1 Absorption apparatus 

The absorption apparatus is shown in Fig. 1. 
It consists of a stainless steel upper section 
containing several gauzes in which the eddies in 
the gas are suppressed, and in which the gas is 
forced through a contracting section, into a flat 
velocity profile. The lower end of this section 
just fits into the top of the cylindrical absorption 
column (glass) in such a way that a slit of 0-4 mm 
width is left; through this slit the liquid enters 
into the wetted wall tower. 

The radial velocity distribution of the gas was 
measured at several distances downstream the 
entrance by means of a hot wire anemometer 
under conditions where the average gas velocity 
and the liquid film surface velocity were equal. 
Figure 2 shows for one example that the velocity 
profile was practically flat, except for a small 
boundary layer effect near the top, which was due 
to wall friction in the contracting section. This 
effect can not, in principle, be entirely eliminated. 
It was possible to fulfill the (arbitrary) requirement 
that the maximum deviation from the average 
gas velocity be reduced to 10 per cent within 
4 of the total height. (See Fig. 2). 

The hot wire anemometer was also used to 
determine the degree of turbulence of the gas 
inside the wetted wall column. From _ these 
measurements it was verified that eddy diffusion 
could be entirely neglected with respect to 
molecular diffusion. 
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Fic. 2. Velocity distribution of the gas for 
v, = 41-4 cm/sec. 


VOL. 2 
ll | 
959/60 
ate 
164 
| 
8 
| 
|| | 
LIQUID OUT ~ 
63 


W. A. Dekker, E. Sworck and H. Kramers 


Whereas the height of the wetted wall column 
(13-6cm) was fixed by the boundary layer 
entrance effect, the inside diameter of the absorp- 
tion column was chosen such that a reasonably 
accurate interpretation of gas absorption results 
from gas analysis data alone would be possible. 
With a diameter of 3-6cm the percentage of 
NH, absorbed by water could be expected to lie 
between 20 and 35, depending on the contact 
time between gas and liquid. 


2.2 Additional equipment 


CO, absorption—For the measurements of the 
absorption rate of CO, into the water film the 
pure CO, gas was saturated with water vapour 
at the temperature of the absorption experiments 
and was then led through the absorber. A very 
small amount of detergent had been added to the 
water in order to prevent surface rippling. The 
CO, absorption rate was found from the decrease 
of the CO, volume at constant pressure which 
was measured by means of the “soap bubble ” 
method [7]. 

NH, absorption -The absorption experiments 
of NH, into water from nitrogen were carried 
out at 25°C. The gas flows of NH, and N, were 
measured separately by means of calibrated rota- 
meters and then mixed and fed to the absorber. 
The incoming and outgoing gas streams 
could be sampled for NH, analysis. The NH, 
content of the outgoing liquid was determined 
by chemical analysis of a number of samples 
during a run. Thus a NH, mass balance could 
be set up around the absorber for every experi- 
ment. 

Also the electrical conductivity of the outgoing 
liquid was recorded in order to verify the steady 
operation during an experimental run. 

NO, absorption—A differential photometer was 
built for the accurate measurement of the NO, 
concentration in the gas mixtures. Its operation 
was based on the selective light absorption of NO,. 
The value of the light absorption coetlicient k, 
for A = 5461 A was determined to be 1-22 em=! 
atm=! by chemical calibration [1]. With this 
photometer the NO, concentration of both inlet 
and outlet gas could be accurately measured at 
frequent intervals. Special precautions had to 
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be taken in order to maintain the temperature 
and the pressure of the gas at constant known 
values because of the sensitivity of the dimeriza- 
tion equilibrium constant of the reaction 
2NO, N,O,. 

During a run liquid samples were taken for 
analysis of HNO, and HNO, in the outgoing 
liquid. In this way, a material balance could 
be set up around the absorber, expressed in 
equivalent moles of NO,. Also in these experi- 
ments the electrical conductivity of the outgoing 
liquid was recorded. 


or THE EXPERIMENTAL 
ResuLts 


3.1 The CO, absorption experiments 


The experimental absorption rates of pure 
CO, gas into the laminar falling water film were 
compared with the rates predicted from the 
“ penetration ” theory for physical absorption. 
According to this theory (see e.g. [7] and [8]), 
the amount of gas absorbed by the liquid with 
initial zero gas concentration will be, per unit 
surface area and after a contact time +: 


d*v, 


3. Discussion 


In a plot of m(r) vs. y/[(h — Ah)/o,] this rela- 
tion would give a straight line through the origin 


with a slope of 2C, y/(D 7). In Fig. 3 such a 
a line has been drawn for 2C,4/(D 7) = 
8-04 kg m*sec’/*, which was found by 
Nusstnc [7] for absorption of CO, at 20°C and a 
CO, pressure of 760mm Hg. The experimental 
points obtained in this investigation under the 
same conditions have also been indicated in 
Fig. 3, and it is seen that the agreement with 
NJSING’s results is very good. 

Also in these experiments it proved necessary 
to reduce the total film height h by the amount 
Ahk which is the height of the reflection wave 
formed above the liquid level of the collecting 
reservoir (1-2. cm in these experiments). In this 
part of the film practically no absorption occurs 
as has been shown in previous work [7, 8]. 

From the CO, absorption measurements it 
appears that the hydrodynamic conditions of the 
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falling liquid film agree well with the theoretical 
model. 


Me, 


Fie. 3. CO, absorption : comparison between 
experimental results and theory. 
A experiment 
theory 


3.2. The NH, absorption experiments 


The flow conditions of the gas were separately 
investigated by absorbing NH, from N, into 
water. Since the solubility of NH, in water is 
high, the absorption rate is mainly determined 
by the diffusion of NH, in the gas phase. If one 
assumes the NH, concentration in the gas at the 
interface, c,,, to be constant the decrease of the 
NH, concentration can be mathematically pre- 
dicted for laminar piston flow of the gas from the 
following formula [9] : 


(2) 


Co — n 


where 8, are the roots of the Bessel function of 
the first kind of order zero. 

In applying this theory to the present case of 
NH, it had to be taken into account that the 
liquid-side resistance is not entirely negligible. 
The value of ¢,, at the interface can be calculated 
for the case where a stagnant gas with concen- 
tration cy and a stagnant liquid are suddenly 
brought into contact with each other. It is then 
found that c,, =f ¢), where f can be regarded as 
the fraction of the total diffusional resistance 
which is due to the liquid phase. In this way, 
f was calculated to be about 0-07 for the system 
N,-HN, — water and under the conditions of 


these experiments. 


Figure 4 shows a graphical representation of 
equation (2) for f = 0 and f = 0-07. When the 
experimental results on the NH, absorption had 
to be plotted in this graph, a difficulty arose as 
to the correct value of the active film height 
and the contact time 7. Since the mass transfer 
is mainly determined by diffusion in the gas 
phase, it seems unlikely that the area of the 
reflection wave just above the liquid level would 
be entirely inactive as was the case with absorp- 
tion of pure CO, into water. Therefore, the 
experimental values were calculated on the basis 


DT 8,)" 
Or/(R-8,)" 


06, 


Or/(R-8,)? x10 


20 30 


Fic. 4. NH, absorption; comparison between experi- 
mental results and theory. 


of the total time of contact (r = h/v,) as well 
as on the reduced time = (h — Ah)/v,). 
Moreover, the accuracy of the results was limited 
by some uncertainty with respect to the mole- 
cular diffusivity of the system NH, — N,*). 

Nevertheless the results plotted in Fig. 4 fit 
the predicted data sufficiently well for the purpose 
of these experiments. It can be concluded that 
transport by turbulence in the gas flow was 
negligible and that the deviation from the flat 
velocity profile at the top of the column had no 
appreciable effect. 


3.3 The NO, absorption experiments 

The absorption rate of NO, into water was 
measured for various contact times (0-2—0-4 sec), 
inlet concentrations of the nitrous gases (3-15 


*A value of 25 x 10-5 m®/see was used for 25°C and 
760 mm Hg. 
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Table 1. 


Main results of the NO, absorption experiments 


Us 


(10-2 m/ 


> 


sec) 


(10-2 atm) 


(m mole min) 


(Peo 


33-8 
33-8 
34:8 
50-3 
641 
O43 


3-16 
2-93 
2-10 
346 
317 
3-42 


3-01 
3-03 
3-02 
2-99 
2-98 
3-09 


340 
35-0 
505 
65-0 


out = 


5-82 
5-96 
5-76 
6-18 
5-69 
5-66 


6-86 
5-74 
6-90 
640 
799 
7:58 


34-5 
34-2 


9-03 


118 
10-9 
12:3 
17 
10-6 
13-4 
13:5 


26-1 
25-2 


per cent) and operating temperatures (25-35°C). 
The main results have been collected in Table 1; 
these lead to the following general conclusions : 


(a) the absorption rate (®,,),, expressed in 
equivalent moles of NO,, is proportional to the 
average partial pressure of N, O, in the bulk of 
the gas (see Fig. 5) ; 

(b) if the inlet concentration of the gas is kept 
constant the absorption rate only slightly depends 
on the contact time ; 


(}m)e(25°C) 
o(¢,), (35ST) 


=e 
32 


Fic. 5. NO, absorption rate as a function of the average 
N,O, partial pressure in the gas. 


(c) the absorption rates are considerably lower 
than would be the case if only the diffusion in 
the gas phase were controlling ; 

(d) the absorption rate, calculated from the 
photometric analysis, is always 30-60 per cent 
higher than that obtained from the chemical 
analysis of the liquid [(®,,),]; 

(e) a change of the operating temperature from 
25°C to 35°C has only a small effect on the absorp- 
tion rate ; 

(f) mist was observed in the gas outlet when 
the eNO, concentration exceeded about 9 per cent. 

From this it can be concluded that N, O, 
rather than NO, plays an important role in the 
absorption process. The deficit in the material 
balance, as mentioned under (d), can be explained 
by the occurrence of a reaction between NO, 
and/or N,O, and water in the gas phase in the 
neighbourhood of wet surfaces. It was estimated 
that at least 93 per cent of this reaction occurred 
in the equipment (tubing and drying column) 
between the gas outlet of the absorber and the 


(°C) (sec) (10°? atm) 
25 0-40 0-29 0-871 
25 0-30 O-882 
35 0-39 0-18 
25 0-27 0-30 0-905 
25 0-21 0-31 0-931 
25 0-21 0-33 0-927 
25 O-s1 O-856 
35 0-39 0-56 O-882 
25 0-27 O-83 0-902 
10 35 0-26 0-62 O-915 
11 25 0-21 O85 0-910 
12 25 0-21 0-81 0-915 VOL 
13 0-39 1-56 0-845 11 
4 0-40 1-52 0-858 1959 
15 343 0-40 9-11 1-61 0-840 
16 a4 0-39 9-15 1-61 0-850 
17 35-6 0-38 9-39 115 0-866 
18 50-8 0-27 9-04 1-63 0-883 
19 65-4 0-21 8-57 1-54 0-899 J 
20 50-5 0-27 15:3 3-16 | 0-864 
21 OT 0-27 15-3 3-30 0-866 
a 
+ + -+—+-] 
2 TTT | 
s + + + + +— + 
6 26 24 28 
66 


The rate of absorption of NO, in water 


photometer, and at the very most 7 per cent in 
the absorber itself. Therefore the photometric 
analysis of the outgoing gas could not be used for 
calculations and the experimental absorption 
rates were computed from the chemical analysis 
of the outgoing liquid. 

The effect, described under (c) suggests that 
the liquid phase also presents a resistance to NO, 
transport. The small sensitivity towards the 
contact time mentioned under (b) might indicate 
that this diffusional resistance is accompanied by 
a rapid reaction between NO, and/or N,O, and 
water. 


4. INTERPRETATION OF THE MEASURED 
AxpsorptTion Rates or NO, into WATER 


For the analysis of the experimental results we 
will use the following model, which is similar to 
that applied by Wenpew and Picrorp [6]: 


(a) Below the entrance of the wetted wall tower 
NO, and N,O,—in continuous equilibrium with 
each other—diffuse from the buik of the gas to 
the gas-liquid interface ; 


(b) N,O, is the only species which dissolves in 
the liquid ; 

(c) the diffusion of N,O, into the liquid is 
accompanied by a rapid pseudo-first-order reac- 
tion between N,O, and water. 


A number of NO, and/or N,O, molecules will 
probably also react with evaporated H,O 
molecules near the gas-liquid interface. How- 
ever, it could be assumed that the influence of the 
gas phase reaction on the NO, absorption rate 
was sufficiently small to be neglected. Moreover, 
complications arising from the liberation of NO 
were assumed to be negligible for the conditions 
of these experiments because of the slow decom- 
position rate of HNQ,. 

Because in particular the physical solubility 
of N,O, in water and the rate constant of the 
reaction mentioned under (c) were not known 
the following calculations were made. 

The diffusion of the NO,—N,O, mixture in 
the gas phase was regarded as the diffusion of 
one fictitious component “ eNO,.” A derivation 
of the diffusivity of this component as function 


of Dyo, and Dy, ,o," Was obtained by assuming a 
steady State diffusion of NO, and N,O, molecules 
through a stagnant gas film of N, (cp. Roperts 
[10}). 

The following differential equation was then 
set up for the diffusion of “eNO,” from the 
cylindrical gas space to the wetted wall: 


= 


dc 
Ty 


(3) 


1 (e/ey) 


(4) 


dd dp? %» 
where p= 1r/R*, = R — and ¢ = Dz/(R*)*v" 


Two of the three boundary conditions are : 


at and 0<p <1, (5) 
and 

c# © at p=0 and d>0. (6) 
At the gas-liquid interface (p = 1) the third 
boundary condition has to account for the above 
hypothesis concerning the N,O, dissolution and 
reaction. The penetration of N,O, molecules in 
the liquid film accompanied by the rapid pseudo 
first order reaction 


N,O, + H,0 “+ HNO, + HNO, 


results in the following equation for the absorption 
rate at the interface : 


(9,,"), = 2C, Vk’ D,) (7) 
As was shown by Danckwerts [11] equation (7) 
is a good approximation as long as k’ z/v, > 4 
which is the case for the greater part of the 
absorbing film in our experimentst. 

At the interface the rate of transport of eNO, 
in the gas phase is given by: 


* Dyo, and Dy,o, were taken to be 1-40 and 0-98 x 10-5 
m? ‘sec respectively at 25°C and 760 mm Hg [2]. 

tStrictly speaking, Danckwerrs’ derivation is valid 
for a constant value of the interfacial concentration in 
the liquid C,. In the present experiments C,, decreases 
somewhat with increasing z but detailed calculations 
have shown here that equation (7) still is applicable with 
good accuracy. 
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The combination of equations (7) and (8) gives 
the third boundary condition : 


Vik D) 


for p=1 and (9) 


dp 


The relation between ¢,, (the eNO, concentration 
in the gas at the interface) and C, (the saturation 
concentration of N,O, in the liquid, which is in 
equilibrium with the gas) follows from Henry's 
law : 


Co = H T (10) 


and the equilibrium between cy, and ¢: 


= 


1/2¢ — K,/8 [4/(1 + 8e/K,) —1]=f(e) (11) 


Substitution of (10) and (11) into (9) finally 
gives: 


(12) 


(¢/€9) AG) 
op Q 
where @ is a dimensionless parameter defined as : 


Table 2. Results of the numerical solution of equation (4) 


T = 


25°C 


3-69 L 
| 


IV 


625 | 


= OO1 


0-935 
0-920 
0-907 
O-895 
O-S84 
O-873 
O-865 
O-853 
OS44 
0-827 


O-o42 
O-929 
0-906 
O-895 
O-S85 
O-S867 
O-858 
O-850 
O-s42 


0-953 
0-932 
0-922 
Ogle 
0-903 
O-SST 


0-500 
0-625 
OSTS 
1-000 
1-125 
1-250 
1375 
1-500 
1-625 


| 
| 


0-046 
O54 
O-901 
O-883 
OST4 
0-866 
O-S858 
O-850 


O-967 
O-958 
0-926 
0-920 
O-907 
0-900 


0-928 
O-S873 
O-s62 
O-S41 
Oss2 


Q = 002 


0-955 
0-933 
0-923 
0-905 
0-806 
O-880 
O-ST2 
0-865 


O-375 
0-500 
0-625 
O-S875 
1-000 
1-125 
1-250 
1-375 
1-500 
1-625 


0-970 
0-954 
O-D47 
0-940 
0-927 
0-921 
0-909 
0-903 


0-936 
O-924 
0-903 
O-893 
O-S87T5 
0-866 
0-858 
0-850 


| 
; VOL 
| | T = 35°C ll 
Series I aa I | II | Ill 
103 119 | 240 | 372 
100 ¢ t/t | | | | 
0-954 
| 
0-932 
| | 0-922 
0-913 
0-904 
| 0-896 
0-888 
0-880 | 
0-873 
| 0-866 | 
| | | 
0-970 0-961 0-980 0-964 
0-962 0-951 0-954 
0-954 0-942 0-968 0-945 
0-947 0-933 0-963 0-937 
0-940 0-924 | 0-957 0-929 
0-932 0-916 0-952 0-921 
0-926 0-908 0-947 0-914 
0-900 | 0-943 0-907 
0-914 0-893 | 0-938 0-900 
0-908 0-886 | 0-933 0-893 
| 0-902 0-879 | | 0-929 0-887 - 
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Since the equilibrium relationship between the 
eNO, concentration in the gas and the N,O, 
concentration in the liquid is non-linear, an an- 
alytical solution of equation (4) with the bound- 
ary conditions (5), (6) and (12) was nearly im- 
possible. Solutions to the differential equation 
were obtained by numerical computation (Table 2) 
for two values of the parameter Q, which were 
chosen by previous trials. 

In order to compare the theoretically computed 
values of @/c, with the experimental results, the 
value of ¢, = D, +/(R*)* had to be determined for 
each run, The contact time 7+ was calculated as 
h/v, and the end effect of the falling film was 
not taken into account for lack of arguments for 
the application of a proper correction. 

The experimental results were plotted as 
é/¢y vs. é, and compared to the theoretical lines 
for Q = 0-01 and Q = 0-02. Figure 6 gives an 


Table 3. 


example of such a comparison for one concentra- 
tion series. 
Then for every run the corresponding value of 


+ 4. 
| 


Fic. 6. Comparison between the eNO, depletion of the 


gas with computed results for Q = 0-01 and 0-02 for one 
concentration series (runs 13, 14, 15, 16, 18 and 19). 


Values of Q obtained from the measurements 


108 
(kmole, m3) 


Uvk D 


102 
| atm sec ) atm 


102 Q 


1-03.10-3 


1-30,10-3 


0-97.10-3 
0-99.10-3 
1-10.10-3 


| 
\ 
Orys 
(25%) 020-01 | 
a 
4 
| } | | 
VOL. 0-80 | 
ll 62466 024 6 6 
959/60 
| 
Run (C) ) 
1 25 1-23 15-8 1 
2 25 oo 15-8 1-4 
4 25 he 10-7 1-2 1-4 | 28 
5 25 = 84 1-6 
6 25 re 84 15 
7 25 2-34 15-1 1-2 
9 25 10-2 1-6 
0-99.10-% 
11 25 8-0 13 
12 25 me 8-0 15 
13 25 3-69 14:5 1 
25 14-7 15 
16 25 14-6 1-3 
18 25 9-9 13 
19 25 7-7 1-2 
20 25 6-25 9-6 10 — 
3 35 1-10 16-9 16 1-6 23 
8 35 2-40 16-2 1-4 
10 35 2-40 10-9 1-6 } sd a6 
17 35 3-72 15-5 1:3 13 27 
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Q was found by linear interpolation between 
Q = 0-01 and 0-02. The results of this interpola- 
tion are shown in Table 3. The Q-values have 
been subdivided into groups, each group represent- 
ing one concentration series and one tempera- 
ture. Reviewing the Q-values of each group given 
in Table 3, it appears that for one temperature 
no significant dependence of the contact time and 
the concentration can be observed. This would be 
expected on the basis of the working theory, since 
Q in equation (13) would only depend on the 
temperature except for a possible influence of the 
concentration on D, and H. 

The change of D, with the concentration is, 
however, relatively small compared to the spread 
in the Q-values, whereas a constant value of H 
was more or less postulated. 

So for each of the two temperatures a gross 
average for Q was deduced from the experiments : 


T = 25°C 
T = 35°C 


Q = (1-3 + 0-2). 10-2 
Q = (1-5 + 0-2). 10-2 


By substituting the values of the known 
quantities* in equation (13) the corresponding 
values of the product H4/k’ were obtained. 
The results are listed in Table 4. 

For purpose of comparison, Table 5 contains a 
number of values for H 4/(k’ D,) obtained from 


* For the diffusivity D, of in water at 25°C a 
value of 1-40 10° m?/see was used [12]. 


Table 4. Final results on the liquid kinetics of 
N,O, absorption 


H y/(k’ Dp 


( 


10-3 kmole =m 


matm — sec m3 atm sect 


11 


absorption experiments performed in different 
equipment. It is seen that the listed values all 
have the same order of magnitude as those found 
in this work. 


5. CONCLUSIONS 


It may be concluded that the experimental 
results of this work are fairly well described by the 
proposed absorption mechanism. Also this theory 
explains that in general the gas phase resistance 
will be predominant at a high concentration of 
the nitrous gases; as the eNO, concentration in 
the gas decreases, the main resistance to mass 
transfer shifts towards the liquid phase. 

In the light of these results it would have been 
more logical to set up a different kind of experi- 
ment in order to study the liquid phase kinetics 
of the NO, absorption and so to obtain more 
accurate values of H y/(k’ D,). As a consequence, 


Table 5. Comparison with results of other investigations 


Authors 


Hy (k D) 
10-3 kmole _™ 


Temperature 
matm sec 


CaupLe and Denpicn* [4] 
CaupLe and Densicn [4] 

and Picrorp [6) 

Wenvet and Picrorp 

SNoECKT 

SNOECK 

SNOECK 


1-10 
250 
0-58 
O54 
0-380 
O77 
O-82 


*Calculated from their experiments for water as the absorbing liquid. 
+These figures result from measurcments by SNoeck in continuation of the work presented in this paper. Here, 


N,O, was absorbed from NO, — N,O, mixtures in a short laminar water jet. 


the work still being in progress. 


‘The data in Table 5 are preliminary only, 
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the NO, absorption experiments have been 
continued with an absorber where a low pressure 
NO,—N,O, mixture is exposed to a laminar 
water jet. Preliminary results of this investigation 
have been included in Table 5. 

In this paper we shall not deal with the calcula- 
tion of the individual values of H and k’. As has 
been shown by Wenpet and Picrorp [6] this 
can be done by combining H 4/k’ values with the 
kinetic studies in connexion with the decomposi- 
tion of HNO, in aqueous solutions, which were 
carried out by Ape. and Scumip [13]. The 


authors think that the treatment by WenpEL 
and Picrorp can be improved by using more 
recent data. This would be worthwhile after 
more accurate values of H y/k' have been obtained 


from new experiments. 
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NOTATION 


¢ = concentration of component to be absorbed in 
the gas phase mole m=? or kg m-% 

C = concentration of gas in the liquid kg m- 
d = internal diameter of the wetted wall tower (dry) 
m 


d+ = 
D = coeflicient of molecular diffusion 


diameter of the cylindrical film surface m 


J = fraction 
h = height of the film m 


= 
H = 
k, = 
kK, 


= gas constant (8-316 x 10) 


- velocity 


height of the end effect m 
Henry coefficient kmole m~* atm-! 
light absorption coeflicient atm-! 


- equilibrium constant of the reaction 


2NO, N,0, kmole m~* 
quantity of gas absorbed per unit area after a 
contact time +r kg m-® 
partial pressure atm 
dimensionless parameter, defined in equation (13) 
cylindrical co-ordinate m 
radius of the wetted wall tower (dry) m 
radius of the cylindrical film = R — dy m 
J kmole~? 
°C °K 

m see} 
surface velocity of the liquid film m sec~! 
co-ordinate in the direction of gas and liquid 
flow m 


temperature 


= thickness of the liquid film m 
= wave length A 
= dimensionless radius r/ R* 


¢= 


®,, = absorption rate per unit area 


contact time (= h v,) sec 
contact time with correction for the end effect 
(= (h — sec 
dimensionless group = Dz, (R*)? v, 

kg m~? 


®,, = absorption rate 
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Advances in Chemical Engineering. Edited by 
T. B. Drew and J. W. Hoorrs, Jr., Academic Press, 
New York, 1958, Vol. II, ix + 338 pp., $9.50. 


HerTeroGeneriry is an inevitable feature of an edition like 
Advances in Chemical Engineering, not only with respect 
to the subjects treated, but also regarding the approach 
by the different authors. Also in this volume it appears 
that the aim to present a “critical summary of recent 
work” (quoted from the editors’ preface) may be 
interpreted rather freely. 

Evidently, the task of the reviewer is to write a number 
of short reviews first and to sum them up, if possible, in 
the end. 

J. W. Wesrwarer finishes his monograph (the first 
part of which appeared in Advances of Chemical Engineer- 
ing, Volume 1) on the Boiling of Liquids (p. 1-3, 33 
references) with chapters on transition boiling, film boiling, 
boiling of subcocled liquids and bumping during boiling. 
Except for the pure film boiling, these phenomena have 
not been successfully explained by theory and the ex- 
perimental material is still very incomplete. This is 
reflected in the text which is clearly written and well 
illustrated, but has to be mostly descriptive. 

Pages 33-80 contain a paper on Automatic Process 
Control (45 references) by E. F. Jounson. Control is 
introduced as a basic concept of chemical engineering and 
control problems are treated analytically, mostly by using 
the frequency response method. This paper has the 
character of an abridged textbook or university course ; 
according to the reviewer's opinion its presentation in 
this volume would only be justified by the fact that the 
application of the servomechanism “language” to 
process control problems is still novel to many chemical 
engineers. 

In Treatment and Disposal of Wastes in Nuclear Chemical 
Technology (p. 81-116, 37 references) by B. Manowrrz, 
a survey is presented of the various methods which can be 
employed for the separation, confinement and disposal of 
radioactive wastes. Clearly, the author has considerable 
experience in the field and so the paper has become a 
useful guide for purposes of orientation. 

G. A. Sorrr and H. C. Weincartner have written 
High Vacuum Technology (p. 117-145, 16 references). 
After a short introduction on the history of the industrial 
application of pressures below 0-1 mm Hg, the authors 
treat the chemical engineering aspects of dilute gases 
(flow, heat and mass transfer, reaction equilibria) with their 
consequences to practical problems. The paper ends with 
a short description of vacuum pumps and gauges. 

Separation by adsorption methods (p. 147-208, 
99 references) by T. VERMEULEN consists of a very valuable 
attempt to unify the theory of various fluid-solids transfer 


operations. Particular attention has been paid to adsorp- 
tion and ion exchange, but application to fixed bed 
absorption or extraction, partition chromatography, 
leaching, drying and regenerative heat transfer are 
possible as well. In the reviewer's opinion, this monograph 
is an important step forward towards the general syste- 
matic treatment of these types of operation. 

Finally, this volume contains thoroughly written 
monograph on the Mixing of Solids (p. 209-324, 73 
references) by S. 5S. Werimpenspaum. The first half is 
devoted to the statistical analysis of the state of mixedness 
and to the treatment of theoretical frequency distributions. 
Chapters on the rate of mixing and the performance of 
various types of dry mixers end this excellent review, 
which is the first of its kind. 

Summarizing, the contributions by Weresrwarer, 
VeRMEULEN and Werpensaum in particular correspond 
to the aim set by the editors, whereas the other papers 
have more the character of handbook or textbook 
articles. The whole collection certainly deserves a place 
in every chemical engineering library. 

The reviewer has been struck by the fact that on the 
total of about 310 literature references 250 are U.S. 
publications ; of the remaining 60, 40 are in the English 
language and about 20 in foreign languages, including a 
reference to Poisuille (1838) and Knudsen (1909), One 
might wonder whether the contributions from the different 
parts of the world to advances in chemical engineering 


are in the same ratio. 
H- Kramers 


Second Symposium on Coal Preparation. Department 
of Mining, University of Leeds, 1958. xiii + 513 pp., 20s. 


Tue Symposium was supported by the Coal Preparation 
Plant Association and held in October 1957 at the Depart- 
ment of Mining of the University of Leeds. Eighteen 
papers were given by authors from Britain, Holland, 
Germany, the U.S.A. and South Africa, and the papers 
are reproduced in the book with the discussion on each 
and summaries by rapporteurs. 

The range of subjects is wide, from the attachment 
between bubbles and particles to the marshalling of railway 
wagons. An idea of the scope of the book may be given 
by dividing the papers into groups according to subject 
and approach. There are several accounts of practice in 
particular fields, including briquetting in the Ruhr, 
coal preparation in Australia, the choice of a material 
for use in suspensions in South African dense medium 
washers, and a general survey of present and future pro- 
blems in Britain. 
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The subjects of the laboratory researches include the 
use of permeability and suction potential methods in a 
study of the retention of water by fine coal, the angle of 
contact between coal and air bubbles, and the attachment 
of bubbles to coal particles in froth flotation. An interest- 
ing discussion of the last paper showed some differences of 
opinion about the relative importance in practice of the 
two collision and gas precipitation—in which 
bubbles become attached to particles, A better under- 
standing of this question should lead to improvements 
in the design of flotation plant. 

Three papers on coal preparation plant deal with jig 
and with dense medium methods, and cover the aspects 
of testing, prediction of performance and representation 
of results. 

Other papers on practical aspects of the subject deal 
with the drying of coal, problems of sedimentation and 
filtration, and experience of interface-active chemicals. 
These last have interesting possibilities as flocculants, in 
spite of their high cost. 

Most of the papers are devoted to the improvement of 
processes that are already well established in the coal 
industry rather than to new developments. The book will 
certainly be of interest to readers familiar with the coal 
industry. Though most of the processes and methods are 
familiar in other connexions, chemical engineers in general 
will no doubt find interest in their application in an 
industry which works on a large scale with processes that 
must be cheap. In general, the field of coal preparation 
offers to the chemical engineer great scope but perhaps 


rather less example. . 
M. J. G, 


Automatic Measurement of Quality in Process Plants. 
Proceedings of a Conference sponsored by the Society of 
Instrument Technology, held at Swansea, Sept. 1957. 
Butterworths, London, 1958. xi + 320 pp., 50s. 


Tue papers contained in this volume form the edited 
proceedings of the 1957 Summer Conference of the Society 
of Instrument Technology and they provide a broad 
survey of the various techniques available for the measure- 
ment and control of quality in process plants. 

In all, twenty-five papers are presented and these are 
divided into six main sections. The first section deals 
with the adaptation of laboratory techniques to plant 
measurement, and a variety of instruments designed to 
perform repeated analytical measurements, such as tit- 
ration and distillation, are described. 

The techniques available for the analysis of gases are 
fully discussed in the second section and calorimetry, 
oxygen analysis and moisture content determination are 
included. 

Continuous measurement of pH and diclectric constant 
are then discussed in the third section. The next two deal 
with spectrometric methods and new techniques for 
fluid stream analysis. Infra-red gas analysis, mass spec- 
troscopy, vapour phase chromatography and the use of the 


property of nuclear magnetic resonance are all treated in 
some detail. 

Various physical properties which can be used for the 
measurement and control of product quality, such as 
viscosity and density, are treated in the last section and 
instruments capable of continuously measuring and 
recording these properties are described. 

Although many of the instruments described are 
capable of providing an accurate assessment of product 
quality it appears that many problems remain to be 
solved before they can be used as the measuring element 
in an automatic control loop. For successful control the 
response of the measuring unit has to be very fast compared 
with the response of the plant and the output from the 
unit should be in a form suitable for use with existing 
control equipment, These conditions are usually met by 
instruments which depend on the measurement of a 
physical property such as temperature, pressure and even 
less familiar properties such as diclectric constant, refrac- 
tive index, viscosity and density, but where chemical 
methods are involved the response of the measuring unit 
is usually comparatively slow and often the output is 
intermittent. As an example we could quote the vapour 
phase chromatograph which does not give a continuous 
reading but a series of discrete records at intervals of 
about 10min. This would require an elaborate control 
system incorporating a memory circuit to store information 
from one record to the next, assuming of course that 
plant is sufliciently stable so that the time interval of 
10 minutes can be tolerated, 

Although chemical methods are slow, the information 
they provide is much more complete than that provided by 
physical methods, especially for multi-component systems 
and in difficult or unusual circumstances their use is often 
unavoidable. 

In summary this volume provides a complete survey of 
the latest techniques available for the measurement of 
product quality and should serve as a very useful reference 
book for process engineers and, to a lesser extent, for 


control engineers. W. L. 


F. S. Martin and G. L. Mites: Chemical Processing 
of Nuclear Fuels. Butterworths Scientific Publications, 
London and Academic Press, New York, 1958. x + 242 pp., 
40s, $7.50. 


Tur book gives a broad picture of the present state of 
knowledge in the nuclear processing field and of expected 


future developments. Particular emphasis is placed on 
the basic chemistry involved. Much useful information 
has been collected and summarized including that released 
at the 1955 Geneva Conference but not from the 1958 one 
nor the 1957 Brussels Symposium. 

The aim of giving as broad a picture as possible and 
underlining the basic principles is an excellent one in 
itself, but inevitably means that less attention is paid to 
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the established processes. The book is commendably 
short but could have been lengthened to accommodate a 
more detailed account of current techniques, particularly 
in aqueous processing. 

It is divided into five sections ; Nuclear Considerations, 
Aqueous Processes, Non-aqueous Processes, Effluent 
Disposal and Fission Product Recovery, and Trends in 
Nuclear Fuel Processing. 

In the first and last of these a good introduction to 
nuclear fuel processing and much useful background 
knowledge is given. A discussion of the economics involved 
shows that fuel processing contributes a substantial part 
to the cost of nuclear power, and yet has no counterpart 
in conventional power production. There is little prospect 
of processing being completely avoided, as a high enough 
burn-up of the expensive fissile fuel cannot be achieved. 

The section on aqueous processes is characterized by a 
detailed discussion of the chemistry involved. There is 
little information provided which would enable a plant 
for separating the constituents of irradiated fuel to be 
designed or its performance to be assessed. Solvent 
extraction has been the mainstay of the U.K.A.E.A.'s 
process for separating uranium, plutonium and _ fission 
products since its conception, and gives every indication 
of remaining so for many years to come. It has world- 
wide application in nuclear fuel processing and hence 
would have justified a more complete discussion. 

Anexcellent survey of the possible non-aqueous processes, 
indicates the paucity of information available. They are 
usually high temperature processes giving partial removal 
of fission products. Although they are simple in concept 
and often leave the fuel in a form suitable for direct 
refabrication by remote handling, they possess technical 
difficulties. A possible application is to the reprocessing 
of fast reactor fuel elements of simple design. 

The problem of radioactive effluent disposal is likely to 
become much more acute in the future. The possible ways 
of dealing with this are well summarized and placed in their 
proper perspective by a discussion of present practice and 
the basic techniques. 

This is a book for the person wishing to acquire a broad 
outline of nuclear fuel processing who is particularly 
interested in the fundamental chemistry involved. It is 
a pleasure to see an English book join the American ones 
on nuclear fuel processing, since substantial contributions 
to this facet of atomic energy have been made by workers 
on both sides of the Atlantic, and anyway, like most 
English books it is more reasonably priced than _ its 


American counterparts. S. 


Scale-up in Practice. Edited by R. Fiemme. Chapman 
and Hall, London, and Reinhold, New York, 1958, 134 pp., 
36s. 


THERE are examples in industry of processes which have 
been successfully taken straight from the test-tube scale 
on the laboratory bench to full size industrial production 


of hundreds of pounds. However, these examples are the 
exception rather than the rule, and more often the basic 
process is developed by a chemist or a physicist at the 
bench and then the chemical engineer steps in and says 
he must build one or even two intermediate scale pilot 
plants to study certain parts of the process before he can 
design full scale constant flow apparatus. This may be 
because the chemist has done a batch process, while the 
production engineer wants to operate a continuous process, 
but it is also very likely to be because of the inherent 
difficulties of scale-up which arise from the changing 
surface to volume ratio, and the fact that different parts 
of the process such as surface and volume chemical reaction, 
heat transfer by natural and forced convection, and 
radiation and turbulence and mixing in fluids have different 
scale-up similarity numbers which cannot all be satisfied. 
Sometimes scale-up even means that none of them can be 
satisfied because the pressure for output on the large plant 
is too great. The commonest difficulty however arises 
from the fact that the chemist says he wants equal time for 
his reaction on the small and the large scale, while the 
physicist wants equal Reynolds number. 


This book is a series of papers presented at a Symposium 
of the American Institute of Chemical Engineers covering 
what might be called a study of the need for pilot plants : 
a not very mathematical treatment of scale-up theory, 
and its limitations based mainly on Damkohler’s treatment, 
but including also Danckwerts’ study of residence time 
distribution which is extremely important wherever 
reaction time comes in. ‘The main conclusion from this 
section is, however, that we do not know enough about 
microscopic mixing which can often be the factor even 
when macroscopic mixing is similar. The next chapter 
deals with the use of analogue computors for solving the 
differential equations of scale-up by simulating them with 
voltage measurement on an electric circuit. A_ final 
answer has not been reached as to where an analogue 
computor is best used and a digital computor, but both of 
them suffer from the disadvantage that they are only 
accurate to the accuracy with which one knows the 
physical properties of the system and can represent its 
geometrical shape in some simple way. However, a good 
case is made for the use of the analogue approach in spite 
of the comparatively high cost. The next paper on pitfalls 
in scale-up gives some very interesting examples of cases 
where scale-up did not go according to plan, particularly 
in the petrochemical industry and petroleum refining, due 
to impurities and accumulation of impurities, and even 
incomplete communication between different departments. 
The next paper emphasizes the need for applying technical 
economics as early as possible in a research and develop- 
ment process and suggests ways of doing this. The sixth 
paper is more concerned with business economics by which 
is meant the factors which help a company to decide 
whether to build facilities for a new product or not. These 
are concerned with the prediction of the size and ground 
of the market, and the probable share of the market, plant 
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location, transport and distribution costs. One might 
consider this chapter a little outside the ordinary use of 
the word scale-up, but it is nevertheless clearly an important 
subject and the more these predictions can be made 
accurately, the better served the industry of a country 
will be. The last chapter deals with organizing for scale-up 
which might also be called administrative arrangements 
to ensure that the chemical industries continue to innovate. 

The book covers a wide field, but obviously none of the 
authors of the chapters have been able to do more than 
give a few practical examples, since it is the record of a 
one-day Symposium. Within this limitation, however, it 
is worth reading by those in administrative posts in the 
chemical industries who are concerned in any way with the 
introduction of new processes and by the more junior 
people who are trying to see that processes promising on 
the laboratory bench have a fair chance of growing to 


commercial status. 
M. W. Turing 


E. Hara, J. Prox, V. Friep and O. Vapour- 
Liquid Equilibrium. (2nd Ed.) 1954. Translated into 
English by G. Standart. Pergamon Press, London, 1958. 
xviii + 402 pp., 90s. 


Tuts is an English translation of a book previously 
published in 1954 in Prague, but in preparing the translation 
the authors have brought their extensive index of published 
vapour- liquid equilibrium data up to date. The book is 
addressed primarily to “ workers in the chemical industry 
who have to deal with problems of distillation and rectifica- 
tion.” It summarizes a very large volume of data and many 
publications that pertain to the thermodynamic properties 
of solutions of non-electrolytes and their implications for 
phase equilibrium between vapour and liquid, and is 
probably the first book of this kind. While other catalogues 
of phase equilibrium data have been available, notably 
the two books by J. C. Cav (1950 and 1956) and his 
co-workers and the collections of azeotropic data by 
Lecat (1918) and Horsiey (1947), no previous publica- 
tions of which this reviewer is aware have brought 
together in a single volume as much of the important 
information derived from thermodynamic principles and 
developments in experimental technique, in addition to 
extensive references to modern sources of experimental 
data. American readers will be especially grateful for the 
many references to data in foreign journals such as 
Chemicke Listy and the Collection of Czechoslovak Chemical 
Communications, including important contributions from 
the Institute of Chemical Technology in Prague, where the 
authors and the translator are staff members. 

The first part of the book contains four chapters devoted 
to general thermodynamic principles, properties of ideal 
and real solutions of non-eclectrolytes and the methods 
that have been found useful for representing equilibrium 
data for real solutions. A second part is composed of five 
chapters on various aspects of laboratory technique used 


in the measurement of vapour—liquid equilibria, such as 
criteria of purity of pure substances ; measurement of 
temperature, pressure, vapour pressure, and boiling point ; 
and direct determination of equilibrium compositions. 
Finally, an alphabetical list of references to over a thousand 
original sources of experimental equilibrium data is 
included, the authors’ literature searches having extended 
to articles published by February, 1957. Many of the 
references probably have not been known to readers 
heretofore. 


The sections on general thermodynamic relations, 
including the Gibbs-Duhem equation and its various uses 
for representing activity coeflicients in non-ideal solutions, 
are based on material that will be very familiar to most 
readers, but it is nevertheless convenient to have such 
material readily available. This is especially true when it 
is associated with a very good summary of the rapidly 
expanding information on empirical equations for excess 
free energies of non-ideal mixtures, such as the polynomial 
equations of Woun, van Laar, Marouies, Repricn- 
Kisrer and others. The authors point out that most of 
the equations in popular use can be derived from Wohl's 
general expressions by various assumptions about empirical 
constants, Several ways for evaluating the constants from 
experimental data are explained in detail and many 
instructive numerical examples are included. However, 
no clear expression of the possible superiority of one or 
another form of the general expansions is made —- probably 
because no such statement is possible at the present stage 
of empirical knowledge. The treatment of methods for 
computing missing compositions from experimental 
observations of temperature-composition or  pressure- 
composition relations for binary mixtures is especially 
well done and clearly reveals the value of such computa- 
tional methods when experimental measurement of 
compositions is difficult or impossible. 

The discussion of activity coefficients of components in 
multicomponent, non-ideal mixtures is also based primarily 
on the formulation of Won, but, unfortunately, the 
authors have employed equations from Wout’s early paper 
(1946) rather than his later one (1953), in which an 
important improvement was made. In Wout’s early work 
the excess free energy for a ternary mixture was represented 
by a polynomial function of the mole fractions which 
typically contained the term ( Ag, + + Age — ©), 
in addition to terms representing contributions from the 
constituent binary systems. The A-constants are derived 
from the binary data but the ternary constant, C, has to 
be found from experimental ternary data if it cannot 
be estimated. The difficulty with this term in the excess 
free energy expansion is that it is ambiguous, for the 
expression is not invariant when the choice of numerical 
subscripts is changed. Thus, if a numerical value of C 
derived from a particular set of data is employed to 
reconstruct the data on the same system, but the subscripts 
1 and 2 are interchanged (and A,, is replaced by Aj,, for 
example), the true values of excess free energy and activity 
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coefficient will not be recovered. The numerical value of 
C in the early Wout equations therefore depends on the 
order in which the pure components are numbered. With 
care in the use of C-values no errors will be made, but 
because of the ambiguity numerical values of C can 
hardly have any physical significance. In his later paper 
Wout replaced the offending term in his equations by 
the alternative expression 7, [}(Ayg + Agy + Ags + 
Agy + Agy + Ayg) — C*), which is obviously invariant 
to all interchanges in the subscripts. Equations derived 
from this modified form are therefore less likely to be 
misused and values of C* can be expected to have some 
basic significance of their own, The same criticism appa- 
rently applies to many multicomponent equations of other 
forms, such as the VAN Laar or SCATCHARD types, and it 
seems regrettable that the present authors have not called 
attention to the difficulties. 

Preparation of a complete summary of vapour-liquid 
equilibrium information of all kinds would be an enormous 
task and would lead to a book at least twice the size of the 
present one. Arbitrary choices of what to include and what 
to omit are therefore inevitable, but it will be a disappoint- 
ment to many chemical engineers that the authors have 


not found it possible to refer more than briefly to the 
technically important literature on K-values and the 
effects of high pressure on equilibria near the critical region. 
The subject is disposed of in two or three pages and not a 
single diagram illustrating critical phenomena and showing 
data for hydrocarbon systems is included. One must 
admit that the illogical division of the subject of vapour— 
liquid equilibrium into one branch dealing with low- 
pressure distillation of organic chemicals (and activity 
coefficients) and another branch dealing with hydro- 
carbons and gases under pressure (and K-factors) has 
become traditional. But a book that summarizes the whole 
subject ought to bring out the unity of the two fields 
rather than this difference. As it is, the information on 
hydrocarbons is included in the list of references to data 
but thermodynamic explanation of many of the associated 
equilibrium phenomena is missing. 

Despite these criticisms the book is an outstanding 
contribution to chemical engineering literature and will 
doubtless be referred to frequently in future studies. One 
hopes that a third edition of the work will appear in a few 
years and will be even larger. 

R. L. Picrorp 
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Annual Meeting of Verfahrens—Ingenieure 


The 1959 meeting of Verfahrens-Ingenicure will take place from 4 to 7 October 1959 at Essen. 


Some twenty papers will be presented and industrial visits will be arranged. 


Further information can be obtained from : 


Die Geschiiftsstelle der VIDI-Fachgruppe Verfahrenstechnik, Frankfurt, M., Rheingau Allee 25, Tel, 770481. 


Second European Symposium on Chemical Reaction Engineering 


On 27, 28 and 29 April 1960 the Second European Symposium on Chemical Reaction Engineering will be held 
in the “ Koninklijk Instituut voor de Tropen” at Amsterdam, Holland. It is being organized by the two Dutch Societies 
“ Koninklijke Nederlandse Chemische Vereniging * and “ Koninklijk Instituut van Ingenieurs” under the auspices of 
the European Federation for Chemical Engineering. 


Secretariat: P. J. Hoftyzer, A M.1.Chem.E., Central Laboratory, Staatsmijnen in Limburg, Geleen (L), Netherlands. 
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